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Abstract. We study the asymptotic behavior of large data radial solutions to 
the focusing Schrodinger equation iwt + Au = — \u\ 2 u in R, 3 , assuming globally 
bounded H 1 (R' i ) norm (i.e. no blowup in the energy space). We show that as 
t — > ±00, these solutions split into the sum of three terms: a radiation term 
that evolves according to the linear Schrodinger equation, a smooth function 
localized near the origin, and an error that goes to zero in the H 1 (R 3 ) norm. 
Furthermore, the smooth function near the origin is either zero (in which case 
one has scattering to a free solution) , or has mass and energy bounded strictly 
away from zero, and obeys an asymptotic Pohozaev identity. These results are 
consistent with the conjecture of soliton resolution. 



1. Introduction 

In this paper we consider the (forward-in-time) asymptotic behavior of global 
energy-class solutions u : [0, +00) x R 3 — > C to the focusing cubic 1 non-linear 
Schrodinger equation 

iu t + Au = F(u) (1) 
in three dimensions, where F(u) := — \u\ 2 u is the focusing cubic non-linearity. 

Non-linear Schrodinger equations such as (1) have been intensively studied and we 
review 2 some of the results here. The equation (1) is locally well-posed in the energy 
class iJ^R 3 ) (see e.g. [10], [45], [41]; in fact it is locally well-posed in H S (H?) for 
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4t is likely that at least some of the results here also extend to other dimensions and exponents. 
However, we need to be in three and higher dimensions in order for the fundamental solution 
e ltA to be "short range" (decaying faster than < _1 ), and the algebraic nature of the cubic non- 
linearity is also somewhat convenient. Also much of the asymptotic analysis is restricted to the 
L 2 -supercritical regime p > 1 + 4 in order for global Strichartz estimates to be useful. 

2 As this is a vast field, we cannot hope to come close to an exhaustive description of results 
here, and the references given are not intended to be complete. We refer the reader to the books 
[10], [45], [6], [41] for a more detailed survey. 
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all s > 1/2 [12]). It also enjoys the conservation of mass 



MASS(u(t)) := 




and conservation of the Hamiltonian 

ENERGY(u(t)) := [ \\Vu{t,x)\ 2 - hu{t,x)\ 4 dx, 
iR 3 2 4 

however these two conservation laws are unfortunately not enough to control the 
i/ 1 (R 3 ) norm; the Gagliardo Nirenberg inequality 

J M 4 <C{J \Vu\^{ J M 2 )^ (2) 

is not sufficient unless the mass 3 is sufficiently small. Indeed, if the Hamiltonian is 
negative then one has blowup in finite time from the virial identity [21], [31] (see 
[10], [45] for further discussion; see also Corollary 1.5 below). 

In this paper we will consider the asymptotic behavior (as 4 t — > +oo) of solutions 
which do not blow up in the energy class; more precisely we consider solutions 
where we assume a priori the uniform H 1 bound 

sup \\u(t)\\ R3 <E (3) 
te[o,+oo) v ; 

for some < E < +oo (note we do not assume any smallness assumption on E; we 
also allow the possibility that solution could blow up at some finite negative time) . 
This condition has appeared in other work on this equation, see in particular [5]. 
In particular (by (2)) we assume that the mass and energy of u is finite. 

Bounded energy solutions will occur for instance when the product of the mass and 
energy of u (which is a scale-invariant quantity) is sufficiently small, thanks to the 
Gagliardo-Nirenberg inequality (2); another option is to modify the non-linearity 
F{u) for large u by replacing it smoothly with a non-linearity which behaves like 
0(|w| p ) for some p < 1 + |, as this allows the Gagliardi-Nirenberg inequality ar- 
gument to derive (3) even when the energy and mass are large (as we shall see, 
it is the low values of u which will dominate our discussion, the high values being 
irrelevant except in order to establish a local existence theory). 

With the a priori assumption (3) we know that u is globally well-posed in if 1 (R 3 ) 
by iterating the local well-posedness argument. But this standard global existence 
result does not reveal very much about the asymptotic behavior of u. For data 
whose energy E is small it is possible to show that the solution u(t) eventually 
scatters to approach a free solution e ltA u + in the energy norm i/ 1 (R 3 ) as t — > +oo, 

3 More precisely, one needs a scale invariant quantity such as (J \u\ 2 )(J |Vu| 2 ) to be small 
at time zero, and then it will be small for all time by combining the above argument with the 
continuity method. 

4 The behavior as t ^ — oo is completely identical, thanks to the time reversal symmetry 
u(t, x) i — ► u(—t, x) of (1), and will not be discussed here. It is however an interesting question 
as to whether the solutions which scatter to a free solution as t — > +oo match at all with the 
solutions which scatter to a free solution as t — > — oo; this is true in the completely integrablc one- 
dimensional situation, as well as for small data, but may well be false for the three-dimensional 
large data case. 
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see e.g. [10], [41] (see also Theorem 1.1 below); note that the potential |u| 2 in F(u) 
is a "short-range" potential (it decays integrably in time if ||u(f)||oo has the expected 
decay of t~ 3 ^ 2 ) and thus does not cause any non-linear corrections (see [31] for more 
discussion). Furthermore, the map u(0) i— » u + is a local homcomorphism in the 
energy space near the origin (see e.g. [41]). 

For large data we do not expect this scattering behavior 5 to the focusing equation 
(1), even assuming the energy bound (3), because of the existence of soliton solu- 
tions (both ground states and excited states) . These soliton solutions are generated 
by finite-energy solutions u(x) to the non- linear eigenfunction equation 

-oju + Au = F(u) (4) 

for some oj > (which is related to, but not quite the same, quantity as the energy E 
in (3)). Such solutions are known to be smooth and exponentially rapidly decreasing 
(sec e.g. [10]), and if one makes the further assumption that u is non- negative and 
spherically symmetric then there is a unique solution to (4) for each u> > [14], 
[41], [3]; we refer to this u as the ground state. There also exist radial solutions 
which change sign, see [2] ; we refer to these as excited states. Note that if u is cither 
a ground state or excited state then u(t, x) :— e lu>t u{x) is a solution to (1); we refer 
to these solutions as the ground soliton and excited soliton solutions respectively. 
These special solutions can be of course modified by the various symmetries of the 
Schrodingcr equation such as scaling, translations, Gallilean transformations, and 
phase rotation to produce solitons of various velocities, energies, widths, etc. 

The analysis of solutions close to a soliton state is by now well understood; the 
ground state soliton for the cubic NLS (1) is unstable due to the existence of 
nearby blowup solutions [1] , [37] , although if one mollifies the growth of the nonlin- 
earity F(u) — \u\ 2 u near infinity enough then one can regain orbital stability [50], 
[51], [11], [22], [23] and asymptotic stability [7], [8], [9], [18], [19] for the ground 
state given appropriate conditions on the nonlincarity; see [10] or [34] for further 
discussion. For perturbations of excited solitons there seems to be far less known; 
one possible conjecture is that generic perturbations of the excited soliton should 
lead either to blowup, or to collapse to a less energetic soliton (or to the vacuum 
state), plus radiation. More recently there has been some results on the analysis of 
perturbations of multisoliton solutions, i.e. superpositions of widely separated and 
receding solitons (see [34], [32]; a similar result for the generalized KdV equation 
is in [26]); while these results do not apply directly to the L 2 -supercritical cubic 
equation (1), they do apply to certain mollified versions of the equation (see [34] 
for further discussion) . We should also mention the work in [46] , [47] on non-linear 



5 One can show however that given any large scattering data £ H 1 (R?), one can find a 
time interval (T„, +00) and an energy class solution u : (T», +00) X R 3 such that limt^+oo \\u(t) — 
e ltA M_|_|| Hl ^3j = 0, basically by solving the Cauchy problem backwards in time from infinity; 
see [10]. But these scattering solutions do not capture all the large data solutions, for instance the 
soliton solutions are clearly not in this class. On the other hand, in the defocusing case every large 
energy (and even certain infinite energy) data scatters to a free solution, while conversely every 
free solution is the asymptotic limit of a nonlinear solution; see [20], [4], [15]. The arguments in 
the defocusing case rely on Morawetz inequalities, which have an unfavorable sign in the focusing 
case; it seems that Morawetz inequalities should still have some value for the focusing case but 
we were unable to obtain any natural application of them in this setting. 
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perturbations of ground and excited states of a Schrodingcr equation with time- 
independent potential. These results however rely either on the variational charac- 
terization of the soliton one is perturbing, or on the fact that the time-dependent 
Hamiltonian —A — \u\ 2 is well approximated by an explicit linear Hamiltonian (e.g. 
the charge transfer model, or the linearization of (1) around a soliton) for which 
the spectral structure is well understood. In the general case we do not have good 
spectral control on this time-dependent Hamiltonian and so these techniques do 
not seem to extend to general large data. 

However, in the special 6 case of the cubic focusing NLS equation iu t +u xx — — 2\u\ 2 u 
in one dimension, which is completely integrable, one can use the methods of inverse 
scattering theory to analyze large data. For instance, if the initial data is smooth 
and rapidly decreasing, then it is known (see e.g. [52], [36]) that the solution 
eventually resolves into a finite number of solitons receding from each other, plus 
some dispersive radiation which goes to zero (although this radiation does not quite 
converge to a linear solution because in one dimension the cubic power is on the 
borderline between long-range and short-range; see [31]). 

Thus one might conjecture that one has similar behavior (resolution into a multisoli- 
ton, plus dispersive radiation, for generic large data (with perhaps some exceptional 
channels such as breathers or slowly divergent solitons for a small unstable set of 
data), or else blowup in if 1 (R 3 )) for other non-linear Schrodinger equations such 
as the ones studied in this paper; we shall refer to this rather informal statement 
as the soliton resolution conjecture. Apart from the results mentioned above when 
the data is close to a soliton or multi-soliton, there are only a few results known to 
support this conjecture. For instance, some progress has been made on controlling 
the growth of higher Sobolev norms for such equations, see [5], [4], [39]. In partic- 
ular for H s solutions to (1) with the bounded energy assumption (3), it is known 
that the H s norms grow at most polynomially in time and that the local H s norm 
stays bounded; these results are consistent with the phenomenon of resolution into 
solitons. The growth of norms between L 2 and H 1 in the infinite defocusing energy 
case have also been studied, see [16], [17], again with polynomial type upper bounds 
on the growth. But to obtain control on asymptotic behavior uniformly in time 
(as opposed to polynomial growth bounds) seems quite difficult. In particular, the 



6 We should caution however that the one-dimensional cubic NLS is L 2 -subcritical, as opposed 
to the three-dimensional NLS which is L 2 -supcrcritical; this seems to play a decisive role in the 
analysis. For L 2 -subcritical equations one has excellent short-time control (in particular, there 
is no blowup), but very poor long-time control even in the defocusing case, whereas for L - 
supercritical focusing equations the situation is reversed. The L 2 critical equation (such as the 
2D cubic NLS) is of course delicate at both short times and long times; see [50], [27], [28], [44] for 
some typical results regarding this equation. 

''This is an oversimplification. In the case when two or more solitons have equal speeds, it is 
possible for "breather" solutions which are periodic or even quasiperiodic, and localized for all 
time to appear. It is even possible to have two solitons of equal speeds and equal heights to recede 
from each other at a logarithmic rate log t in one dimension, although this latter phenomenon may 
be also due to the slow decay of the fundamental solution in one dimension and we do not know if 
it occurs in three dimensions. Nevertheless the generic behavior is that of resolution into receding 
solitons of different speeds, plus radiation, and even in the exceptional cases the non-radiation part 
of the solution is still very well localized and smooth, and has asymptotically constant velocity. 
We thank Jim Colliandcr for pointing out this subtlety. 
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soliton resolution conjecture for non-integrable equations and generic large data 
seems well out of reach of current techniques. 

The case of general large data seems too difficult to analyze at present, so we 
now restrict our attention to the much simpler spherically symmetric case, where 
u(0,x) = u(0, \x\) (and hence u(t, x) — u(t, \x\), by the rotational symmetry of (1) 
and the uniqueness theory). In this case the radial Sobolev inequality 8 

III x MIl~(R 3 ) ~ IMIffi(R 3 ) ( 5 ) 

combined with (3) tells us that u(t) decays away from the origin 9 for all times t, 
and thus we only expect at most a single stationary spherically symmetric soliton 
at the origin, and no other soliton behavior anywhere else. 



Unfortunately even in the spherically symmetric case we were still unable to settle 
the soliton resolution conjecture. However, we have a number of partial results 
which may be of interest. We begin with two rather standard results (perhaps 
already known, though we could not find them explicitly in the literature). First, 
we show that if the solution ever has sufficiently small mass near the origin at a 
sufficiently late time, then one in fact has scattering and no soliton will ever form: 

Theorem 1.1. There exist constants e, R > depending only on the energy E such 
that the following statement is true: if u is any spherically symmetric solution to 
(1) obeying (3), and 

lim inf / \u(t, x)\ 2 dx < e, 

then there exists a function u + G i7 1 (R 3 ) such that 

KmJu(t)-e itA u + (t)\\ Hl{R3) =0. 

Here 5(0, R) is the ball B(0, R) := {x G R 3 : \x\ < R}. 



We prove this easy result (a variant of the small data scattering theory, combined 
with (5)) in Section 4. 

From Theorem 1.1, we see that in order to have non-trivial asymptotic behavior 
(such as the presence of solitons) , we must permanently station some mass near the 
spatial origin as t — > +oo. It is then natural to ask what the asymptotic behavior of 
this mass is. One can show without too much difficulty, even without any assump- 
tion of spherical symmetry, that the solution decouples into a free solution (which 
thus decays near the origin, thanks to local smoothing estimates) and a remainder 



8 This inequality can be easily deduced by applying the one-dimensional Sobolev inequal- 
ity IMI_l=o(PI) < IMIifi (R) to * ne function u(|a;|) := |a;|u(|a;|) and using Hardy's inequality 

IN/MII^rS) <l|v«|| i2 ( R3) . 

9 This explains our previous remark that the large values of u are not particularly relevant for 
our purposes; indeed, we see from (5) that they only make a difference near the spatial origin. 
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which is "weakly bound" in the sense that it is asymptotically orthogonal 10 to all 
free solutions: 

Theorem 1.2. Let u be any solution to (1) obeying (3) (not necessarily spherically 
symmetric). Then there exists a unique decomposition 11 

u(t) = u wb (t) + e UA u+ (6) 

into the weakly bound component u wb and the linear radiation component e ttA u + , 
obeying the energy bounds 

H u +Hffi(R 3 )' SU P H u w&WHffi(R 3 ) ~ 1 ( 7 ) 

v ; te[o,+oc) ' 

and the asymptotic orthogonality conditions 

t KmJu wb ,e itA f) L 2=0 (8) 

for all f G i? _1 (R 3 ), or equivalently that 

t ^J u wb^ eitA f)^=0 (9) 

for all f G H^R 3 ), where 

(/, 9)m ■= if, q)l> + (V/, Vg) L , = f fg + V/ • V5 dx 

is the inner product for the Hilbert space i/ 1 (R 3 ). Furthermore, u wb is a (local-in- 
time) approximate solution to (1) in the sense that 

T^+oo ll( * 9t + A)U ™t - F ( U wb)\\Lim A lT,T+T]xIl 3 ) = ( 10 ) 

for all t > 0. Finally, we have the mass decoupling identity 

MASS(u+) + r lim MASS(u wb (T)) = MASS(u(0)) (11) 

and the energy decoupling identity 

\ J R3 |V«+| 2 + T ^n oo ENERGY(u wb {T)) = ENERGY(u(0)); (12) 

in particular the above two limits exist. Finally, if u is spherically symmetric then 
u wb and u + are also spherically symmetric. 



We prove this theorem (which is another variant of the small data scattering theory) 
in Section 5. We remark that (11), when combined with Theorem 1.1, shows that in 
the spherically symmetric case the asymptotically bound mass limT^+oc M AS S (u w ^(T)) 

10 More precisely, this is an asymptotically bound state at t = +oo. There is also of course a 
notion of an asymptotically bound state at t = — oo, and it is an interesting question (which we 
were unable to address) as to whether these two notions are at all related. 

^This is perhaps not the optimal decomposition of the solution u, as it only isolates the linear 
evolution component of the radiation term, and places everything else (including non-linear self- 
interactions of the radiation term) into the "bound state" term. However, modulo errors which 
go to zero in the energy norm as t — > +oo, this seems to be the correct splitting into radiating 
and non-radiating states; this is confirmed in the radial case by Theorem 1.3 below. 
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is either equal to zero, or is bounded from below 12 by some e — e{E) > 0. As for 
the asymptotically bound energy, see Corollary 1.5 below. 

Theorem 1.2 is consistent with the idea of global solutions to (1) decoupling into a 
bound state plus radiation, but the control on the weakly bound state u w ^ given 
by the above Theorem is very unsatisfactory; for instance, we do not know whether 
u wb a PP roac hes a global solution to (1) in any reasonable sense (the above result 
only shows that for any fixed length of time r > 0, the solution u w ^ approaches a 
solution to (1) on intervals [T, T + r] for T sufficiently large depending on r). In 
particular we do not know whether u ^(t) converges asymptotically to a soliton, 
even if we restrict the time t to a subsequence of times. However, in the spherically 
symmetric case we can prove that u w ^ becomes asymptotically smooth and has 
some decay: 

Theorem 1.3. Let u be any spherically symmetric global solution to (1) obeying 
the energy bound (3). Then for any J > 1 and S > 0, we can decompose the weakly 
bound component u w ^ of (6) further as 

u wb(t) = u b^) +°ffi(R 3 )( 1 ) ( 13 ) 

for allt > 0, where °jji/p£ 3 )(l) * s a time- dependent error whose ij^R 3 ) norm goes 
to zero as t —* oo, and u^{t) is a spherically symmetric function (depending on J, 
5) obeying the symbol-type estimates 

\Viu b (t,x)\<(x)-i-j +s (14) 

for all x G R 3 , t >T, and < j < J, where (x) := (1 + l^l 2 ) 1 / 2 and the implicit 
constants can depend on the exponents 5, J. Also, we have 

sup IHWII H i ( R 3 ) £ 1- (15) 



This theorem is considerably more technical than the previous two, and is the main 
result of this paper. It is proven in Sections 6, 7, 8. As mentioned earlier, the 
assumption (3) can be dropped if we mollify the non-linearity smoothly at infinity 
to decay like 0(|u| p ) for p < 1 + |, although this requires a number of technical 
changes to the argument (e.g. the introduction of the fractional chain rule) which 
we omit. 

Note that the error term in (13) is only controlled in the homogeneous space 13 
_ff 1 (R 3 ) instead of the more natural energy space i/ 1 (R 3 ); our arguments have the 
defect of not being able to control the very low frequency portion of the solution u in 
a satisfactory manner (in particular, we could not prevent the very low frequencies 
of Ufa from absorbing a non-trivial amount of mass in the asymptotic limit) . This is 
also related to why the bounds in (14) only give decay of 0((x) _3 / 2+ ) instead of the 

12 Strictly speaking, this argument also requires Lemma 3.3 to dispose of the possible effect of 
the radiation term e rfA u_|_. The above remark can also be derived in the non-sphcrically symmetric 
case by a standard perturbation analysis of the free solution e ltA M_(- at large times, similar to the 
construction of wave operators for this equation; we omit the details. 

■^Of course, since H 1 controls L 6 by Sobolev embedding, we do get convergence in H 1 (R?) 
on any fixed compact set in space, though again this would be more useful if we could improve 
the bounds on in (14) from 0({x)- :i / 2 +) to the square-integrable 0{{x)- z / 2 -). 
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exponential decay which is to be expected by comparison with the soliton solutions. 
If one could improve the error control to H , and also improve the bounds in (14) 
from 0((x)~ 3 / 2+ ) to 0({x)~ 3 / 2 ~) or better, then the above theorem would say that 
the evolution of (1) in the bounded energy case (3) has a compact attractor in H 1 
(once one projects out the linear radiation e ttA u + ), namely the set generated by 
(14). In principle this makes the evolution of almost periodic, and reduces the 
asymptotic study of (1) to the dynamical systems problem of understanding the 
flow on or near this compact attractor. It also raises the possibility that one could 
analyze the asymptotic behavior of the bound portion of the solution by passing 
to a weak limit at infinity. We caution however that the above theorem does not 
establish that wj-, is in any way a global solution to (1); the best we could do was 
choose itfo so that it obeyed bounds of the form (10) (see also Lemma 10.1 below). 

Thus Theorem 1.3 is still far from satisfactory in a number of ways, though it does 
represent some progress towards the soliton resolution conjecture. We speculate 
that the use of such identities as the virial identity, pseudoconformal conservation 
law, or the Morawetz inequalities (none of which are used in the proof of Theorem 
1.3, which relics almost entirely on an analysis of the high frequencies only) may be 
useful tools to analyze this problem further; we give one example of this in Theorem 
1.4 below, when we use virial-type identities to prove an asymptotic Pohozaev 
identity. 

To put Theorem 1.3 in some context, we may compare it with the elliptic analogue, 
which is as follows. Consider a solution u{x) to the nonlinear eigenfunction equation 
(4) for some ui > 0, where u is assumed to be spherically symmetric and finite 
energy, but not assumed to be non-negative (since that would force u to be the 
unique ground state [14]). To put it another way, we assume that u is a critical 
point of the Lagrangian associated to (4) but is not necessarily the minimizer. The 
analogue of Theorem 1.3 would be then that u is smooth and rapidly decreasing. 
In the elliptic case, both facts follow easily from rewriting the equation (4) as 
u = (u — A) _1 (|u| 2 u) and then iterating this equation to bootstrap the finite 
energy and spherical symmetry assumptions to obtain arbitrary smoothness and 
decay; see e.g. [10]. 

In the dispersive case, where we have (1) instead of (4), we do not have the luxury of 
a smoothing operator such as (uj — A) -1 . However, we do have dispersion, which in 
principle achieves a similar effect (locally in space, at least) given sufficient amounts 
of time. The heuristic justification for the asymptotic regularity of u is that any 
high frequency component of u near the origin should quickly radiate to spatial 
infinity as t — > ±oo, but because u(t) has bounded energy as t — > +oo this radiation 
should eventually decay in time on each fixed compact region of space (cf. the local 
smoothing effect [38], [48] for the Schrodinger equation). Thus if we restrict to 
t > T for T sufficiently large, the total energy of high frequencies that pass near 
the origin for times t > T should be small. The main difficulty here is to ensure that 
the high frequency radiation does not interact significantly with the "bound states" 
of the time-dependent Hamiltonian — A — |u| 2 . This gets easier to accomplish if one 
already has some preliminary regularity result on u, which suggests that we can 
prove the regularity part of Theorem 1.3 by an iteration argument starting from 
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(3), similar to how one proceeds in the elliptic case. It also transpires that the 
further one gets away from the origin, the larger the range of frequencies one can 
classify as "high frequency" (roughly speaking, any frequency higher than (x)~ 1+s 
can escape to spatial infinity either as t — > +00 or t — > —00, depending on whether it 
is outgoing or incoming with respect to the origin) . This explains the decay factors 
in (14). Unfortunately we do not get infinite decay, because we could not find a 
dispersive analogue 14 of the elliptic fact that the smoothing operator (uj — A) -1 is 
highly local. 

As indicated above, the methods used to prove Theorem 1.3 are fairly elemen- 
tary, relying principally on the Duhamcl formula, the Strichartz estimates and the 
dispersive inequality. There is one somewhat exotic ingredient however, which is 
a microlocal decomposition of free spherically symmetric solutions into incoming 
waves (waves moving towards the origin) and outgoing waves (waves moving away 
from the origin). This decomposition is useful because we have a freedom when 
applying the Duhamel formula to evolve either forwards in time (controlling the 
present from the past) or backwards in time (controlling the present from the fu- 
ture, taking advantage of the bounded energy hypothesis (3) on the future). To 
analyze outgoing waves it is easier if one evolves forwards in time (as the wave be- 
comes incerasingly dispersed), but for incoming waves it is easier to go backwards 
in time (in order to avoid the wave concentrating at origin, where the non-linear 
effects are strongest). To exploit this decomposition we will rely heavily on duality 
(and in particular of the unitarity of the free evolution operators e ltA ). This decom- 
position may be of independent interest (for instance, a similar decomposition is 
utilized in [35] to establish global-in-time local smoothing and Strichartz estimates 
on compact non-trapping perturbations of Euclidean space). 

One of the main drawbacks of Theorem 1.3 is the poor control of the decay of u^, 
which is only shown to be 0((x)~ 3 / 2+ ). In particular, we were not able to place 
xu-fo in L 2 , which seems a natural objective as one may then apply the virial identity 
or pseudoconformal conservation law directly to obtain further information on u^. 
We were however able to exploit the weaker decay in Theorem 1.3, combined with 
a modified virial identity, to prove the following statement: 

Theorem 1.4 (Asymptotic Pohozaev identity). Let u is any spherically symmetric 
solution to (1) obeying (3), and let be as in Theorem 1.3. Then we have 

1 r T+T r 

lim (lim sup \- / 4|Vud 2 - 3|u fo | 4 dxdt\) = 0. 



The difficulty is that, unlike the elliptic case, the portion of u near the origin can eventually 
influence even the very low frequencies u far away from the origin by Duhamcl's formula (27). To 
control this influence it seems that we need to control the spacetime integral of the non-linearity 
\u\ 2 u on large regions of spacetime. Using the equation (1) one can get some reasonable control on 
such spacetime integrals, but not infinite decay. Interestingly, in the case of a soliton u = Qe~ lult 
one can get infinite decay on the spacetime integral of \u\ 2 u thanks to the time oscillation; if such 
a phenomenon held for general solutions then we would also get infinite decay of u^. 
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As the name of the Theorem suggests, this theorem is the dispersive analogue of 
the Pohozaev identity 



L 



4|Vw| 2 - 3|w,| 4 dx = 

R 3 



for solutions to the non- linear eigenfunction equation (4), which is proven by mul- 
tiplying (4) against x ■ Vu and then integrating by parts. The above theorem is 
a (rather weak) generalization to the dispersive case, saying that the Pohozaev 
identity is asymptotically verified in an average sense for the bound portion of the 
solution (a similar estimate holds for u w ^, by (13), Sobolev embedding, and the fact 
that both Ufo and are bounded in H 1 ). The standard dispersive generalization 
of the Pohozaev identity is Glassey's virial identity 15 

d t j 2x- Im(Vmt) dx = j A\Vu\ 2 - 3|u| 4 dx, (16) 

which, unsurprisingly, is also obtained by multiplying (1) by x ■ Vu and then inte- 
grating by parts. Unfortunately the finite energy assumption on u is not enough to 
ensure that the integral on the left-hand side of (16) is finite (even with the spher- 
ical symmetry assumption), however if we replace u by and we add a slight 
damping weight of (x)~ 5 then (14) will suffice to make the integral finite. We carry 
out this (standard) procedure and prove Theorem 1.4 rigorously in Section 10, after 
a review of virial identities in Section 9. 

Combining Theorem 1.4 with Theorem 1.1 yields the following lower bound on the 
asymptotically bound energy: 

Corollary 1.5. Let u is any spherically symmetric solution to (1) obeying (3), and 
such that the asymptotically bound mass lim t ^ +oc MASS{u w ^) is non-zero. Then 
there exists a = <r{E) > such that the asymptotically bound energy is also bounded 
from zero: 



lim ENERGY(u h ) = lim ENERGY(u v , h ) 

^+oo u t— >+oo wu 



> (T > 0. 



Proof Let u be as above. By Theorem 1.1, there exists so > and R > 
(depending only on the energy E) such that 

lim inf / \u(t, x)\ 2 dx > e - 
t ^+°° Jb(o,r) 

By (6), (13), and Lemma 3.3 below, we thus have 

lim inf / \uu (t, x)\ 2 dx > en. 
In particular, from Holder we have 

lim^ inf J |u^(t, x)\ 4 dx > c 



15 The left hand side is also equal to dtt f M 2 M 2 dx, however in our applications we have 
nowhere near enough decay to make sense of this quantity and so we will rely instead on the 
left-hand side as stated in (16). 
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for some c = c(so, R) > 0. As a consequence, we have 

lim t vrf^J 3 4 I V %W| 2 - 3\u h (t)\ 4 dx < lim t inf^ 8ENERGY(u h ) - c. 

By Theorem 1.4 we thus obtain 

lim inf 8ENERGY(u h ) - c > 0. 

Observe from (7), (15) that the error °H-i(p_ 3 )(l) m (13) is bounded in i/ 1 (R 3 ) 
norm, and in particular by (2) will be decaying in L 4 (R 3 ) norm. From this it is 
easy to see that 

t lim ENERGY(u h ) - ENERGY(u wh ) = 0, 

and the claim then follows from the above estimates (and (12), to establish existence 
of the limits). ■ 



Corollary 1.5 is a variant of the results in [21], [30], which state that initial data 
(which is cither localized or spherically symmetric) with negative energy will lead 
to blowup in finite time. It is also consistent with soliton resolution, since solitons 
have strictly positive energy in the L 2 supercritical case (see e.g. [10]). 

Remark. Except for Theorem 1.2, our results are restricted to the spherically 
symmetric case (although it seems Theorem 1.1 generalizes fairly easily, see the 
remarks at the end of Section 4), mainly because (5) gives us excellent control on 
where the solution is large (at the origin) and small (everywhere else) . The general 
case is of course far more difficult, but one could speculate that similar results ob- 
tain for general bounded energy solutions to (1). At any given time there should be 
a bounded number of "points of concentration" , and the weakly bound component 
u w ^ should be asymptotically smooth and concentrated near these points, and de- 
cay away from these points. Furthermore these points should move at bounded 
speeds, though it is not clear to the author whether these speeds will eventually 
stabilize to be asymptotically constant. If these points of concentration end up 
receding from each other, then the behavior of the bound state near each concen- 
tration point should behave much like the analysis given here for the spherically 
symmetric case (after a Gallilean transform), in analogy with the analysis in [34]. 
Thus ideally one should be able to reduce the problem of soliton resolution to the 
case where the only concentration point is the origin, though as with the results 
here, to proceed further seems to require some analysis of the dynamical system 
of the compact attractor associated to this concentration point, which the author 
docs not know how to pursue. 



2. Notation 



We use the notation X < Y or X = 0(Y) to denote the estimate X < CY, where 
C > is a constant that can depend on exponents (such as regularity exponents 
j or a), as well as the energy E in (3), but not on other parameters such as time 
parameters t or functions u. 
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We shall abuse notation and write /(|a;|) = f(x) for spherically symmetric functions 
/, i.e. we think of such functions as living on [0, +00) as well as living on R 3 . 

We fix 77 to denote a smooth spherically symmetric function on R 3 adapted to the 
ball 5(0,1) which equals 1 on 5(0,1/2). For any R > use t]r to denote the 
rescaling rjn{x) := r)(x/R) of x, thus vr is adapted to B(0,R) and equals 1 on 
B(0, R/2). We shall frequently use differences such as wr — r\ r or 1 — 7]r to localize 
space smoothly to various annular regions. 

We define the Fourier transform /(£) on R 3 by 

We then define the Littlewood-Paley operators Pj for all integers j by 

(17) 

We record Bernstein's inequality 

II/IIl. ( r», S 2(f - f)j 'll/ll LP(R 3 ) 

whenever j 6 Z, 1 < p < g < 00, and the Fourier transform of / is supported on 
the ball 5(0, CV~) (e.g. / could equal Pjg for some g). This can be proven by 
applying Littlewood-Paley type operators and using Young's inequality, as well as 
the standard bounds for the convolution kernel of the Littlewood-Paley operators. 

We also define the ordinary Sobolev spaces via the norm 

ll/ll*.(R») : = IKO'/II^R., 

and the homogeneous spaces H s (R 3 ) by the norm 

H/lliHR 3 ) : = \Wt\'fW»&r 

We shall unify these two spaces by defining the scaled Sobolev space £/^(R 3 ) for 
any spatial scale R > by 

H/Hffg(R 3 ) : = (/ R3 (ie| 2 + ^) Q |/«)| 2 ) 1/2 

Thus Hf (R 3 ) is the inhomogeneous space H a (R 3 ) , while iP(R 3 ) is m some sense 
the limit of H^(H 3 ) as R — > 00. There is also the easily verified scaling relationship 

ll/«ll^(R 3 )=^ 3/2 " a H/ll^(R3) 

for all R, where Jr{x) := f(x/R). We remark that i?^ is the dual of H^ a , and 
that we have the estimates 

II/H^(R 3 ) £ II/H^(R 3 ) £ H/IIh-(R») £ ^ll/llngcR 3 ) 
when -R > 1 and a > 0, and dually that 

II/H^(R 3 ) * II/H^(R 3 ) * II/H^(R 3 ) £ ^H/H^R-) 
when i? > 1 and a < 0. 

We also have the following handy stability lemma with respect to cutoffs. 
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Lemma 2.1. Suppose that a(x) is a function on R 3 obeying the estimates 

sup |V J a(x)| < AR- j 

for some A > and R > ; and a// < j < J. Then we have the estimate 

/or all — J < a < J, where the implicit constants may depend on the exponents a 
and J. 



Proof We may rescale A = 1 and R = 1. By duality we may take a > 0; by 
interpolation wc may assume that < a < J is an integer. But then the claim is 
clear from the Leibnitz rule and Holder's inequality. ■ 



We shall also need the standard Sobolev spaces 

k 
3=0 

and 

ll/ll^ P (R 3 ) := H^IL^R 3 ) 
for 1 < p < oo and non-negative integers k. 



3. Preliminary estimates 

Fix u solving (1) and obeying (3); to begin with we do not assume spherical sym- 
metry. 

Let e ltA be the propagator for the free Schrodinger equation iu t + Au = 0; this is 
of course a unitary operator on L 2 (and indeed on every Sobolev space H S (R 3 )). 
From the explicit formula 16 

e UA f(-)--^j Il3 ^ 2/4t f(y)dy (18) 

we obtain the standard dispersive inequality 

II^/IUr^^II/ILhr 3 ) ( 19 > 



We ignore the issue as to which branch of the square root of t to take when t is negative by 
hiding this issue under the undisclosed constant C. 
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and as a consequence the global Strichartz inequality (see [24] and the references 
therein) 17 

H^ A /llL ?iS( RxR 3 ) + H^/IL-^RxR-) S Wf\\ LH Wy ( 2 o) 

and hence (since e ttA commutes with derivatives) 

ll el * A /IL?^' 6 (RxR 3 ) + \\ eltA fhrm(RxR, 3 ) ~ H/Hff 1( R 3 )- ( 21 ) 
In particular we observe from Sobolev embedding that 

H e ** A /IL ? L=o ( RxR 3 ) + H e ** A /llLfi«(RxR 8 ) + \\ eUA fh?Ll(IlxIl 3 ) ~ H/llifi(R 3 )- 

(22) 

As observed in [24], there are inhomogeneous versions of these estimates in which 
one introduces a forcing term in a dual Strichartz space. For instance, we have 

ll U IL~,fJi([to,+oc)xR 3 ) + H M llL*L6([t ,+oc)xR 3 ) 

< IK^U^RS) + ll*«t + Au|| 1 . ?w .1.6/. ([toi+oo)xR 3 ) (23) 

for any (Schwartz) function u and any time to- Of course, many other inhomo- 
geneous Strichartz estimates are available, but we isolate this particular one as it 
shall arise specifically in our arguments 18 . 

We also need the following variant of above estimates: 
Lemma 3.1. We have 

H e " A /IL^(RxR 3 ) S H/H^RV ( 24 ) 

Proof This almost follows by interpolating (21) with itself, but the presence of 
the means that endpoint Sobolev embedding cannot be applied directly, and 
we must proceed with a bit more care, using interpolation theory instead. We 
shall use the Littlewood-Paley operators P, defined earlier. Observe from (20) and 
Bernstein's inequality that 

ll^/IL^R^ < 2^||P J /|| L2(R3) ~ 2-'/ a ||P J /|| Al(R , )I (25) 
for any j, and similarly that 

ll^/IL-^^R.) < 2^/ 2 ||P,/|| L2(R 3 ) ~ 2^||P j /|| Al(R 3 ) . (26) 
By Marcinkeiwicz interpolation we thus obtain the Lorentz space estimate 



^ A ^/ll^ lL oo (RxR 3 ) < ||iVI|* 1(R * 



17 In this radial setting there are in fact some extra smoothing estimates available, see [49], 
although for our argument we will not need any such refinements of Strichartz inequalities (nor 
will we need more complicated versions involving for instance X s ' b spaces. The main difficulty 
in this work is not so much the recovery of local regularity, which is fairly easy since we are in 
a sub-critical setting with no derivatives in the nonlinearity, but rather in obtaining sufficient 
long-term control (e.g. decay) of various components of the solution in time, and the relatively 
basic Strichartz and dispersive estimates already seem adequate for this task, at least in the 
three-dimensional case. 

18 Note that these Strichartz estimates are not particularly sharp; we concede a number of 
derivatives and do not always use the endpoint exponents. This is a reflection of the fact that the 
non-linearity is not at a critical power such as the L 2 -critical power p = 1 + ^ or the iP-critical 
power p = 5. 
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and thus by the triangle inequality we have the Besov space estimate 

ll eJtA /IL-^(RxR 3 )^Ell^/H ffl (R 3 )- 
jeZ 

On the other hand, from (25) and (26) we have the estimate 

^||e^ A P j /|| L „ (R 3 ) rft<min(2-^|^|V2 )2 i/2| S | ) ||p. / ||^ (R ^ 

for any measurable set E. Summing this in j, we obtain 

/ H^/II^R-) dt Z |£| 3/4 sup 11^/11^) 
J E J 

and hence (since E was arbitrary) the weak-type estimate 

II^/IL-l-cr^^^^pII^/Uhmr 3 )- 

Interpolating this with the previous Besov space estimate we obtain the result. ■ 

We can extend these estimates from the linear evolution to the non-linear evolution 
(1), provided we localize in time. The key tool here is the well-known Duhamcl 
formula 

rti 

u(ti) = e^-^uito) - i \ e l{tl - t)A F(u{t)) dt (27) 

for solutions to (1) and all times to,h 6 [0, +oo), where we adopt the convention 
that J 4 1 — — f t ° if t\ < to- Of course, it is the time integral in (27) which poses the 
most problems, especially as \t\ — to I gets large. Generally speaking, our strategy 
in this paper will be to apply Strichartz estimates or local smoothing estimates in 
the short term (when \t\ — t\ is small) and to use the dispersive inequality (19) in 
the long term (when \ti — t\ is large). 

We now give the (standard) non-linear local-in-time analogue of the Strichartz 
estimates. 

Lemma 3.2. Let u be a solution to (1) obeying (3). Then we have 



for all 2 < p < 6, and 



for all T > and t > 



for allT>0 and t > 



sup K*)|| iP(R 3 <1 (28) 
te[o,+oo) 



II«IIl ? ^^([t,t + .]xR 3 )^M 1/2 ( 29 ) 
In particular we have 

ll^l^dT^xR 3 )^^) 172 ( 3 °) 
Similarly, we have 

ll^ll^dT^xR 3 )^^) 174 - (31) 
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Proof The bound (28) follows immediately from (3) and Sobolev. To prove (29), 
(30), (31), it suffices to do so for r > sufficiently small depending on E, since the 
claim for larger time spans r then follows by decomposition of the time interval 
into smaller pieces. We now fix T, r, and let X denote the quantity 

X := H U IL?wi- 6 ([T,T+r]xR 3 ) + IHIl?L~([T,T+t]xR 3 ) + W U W L% L~([T,T+r] xR 3 ) ' 

By (21), (22), (24), (27) and Minkowski's inequality we have 

X Z \H T )W^(R') + jl +T ll^(«(*))ll H i(R») dt - 

By (3) the first term on the right-hand side is O(l). For the second term we use 
Leibnitz and Holder to bound 

r-T+T 

J t ||FWt))|| Hl(R3) dt < ^llull^j^^^RajllulLoo^j^+^^S) 

and so by (3) and definition of X we have 

X < 1 + t 1/2 X 3 

and the claim follows by standard continuity arguments (e.g. letting r increase 
continuously from 0) if r is sufficiently small. ■ 

Finally, we record a Ricmann-Lcbesgue type lemma for the free Schrodingcr prop- 
agator e 2tA , which while preserving the H S (H 3 ) norms, sends many other norms to 
zero: 

Lemma 3.3. For any u e iJ x (R 3 ) we have 

lim f (x)- E (\e ttA u a (x)\ 2 + \Ve ttA u (x)\ 2 ) dx = 
for all e > 0, and similarly 

lim [ \e ttA u {x)\ p dx = 
*^+°° JR 3 

for all2<p<6. 

Proof Observe that the expressions in the limits are certainly bounded by some 
quantity depending only on the H 1 (R 3 ) norm of u , thanks to Sobolev embedding 
and the fact that e %tA preserves the H 1 norm. Thus by the usual limiting argument 
(and the linearity of e ltA ) it suffices to verify this lemma for test functions u , which 
are dense in H 1 (R 3 ). But in that case we see from (19) that \e ltA u (x)\ decays like 
0{t~ 3 / 2 ) (with the implicit constants here depending on uq). Since |e l * A uo(a;)| is 
also bounded in L 2 , the claim follows. ■ 



4. Proof of Theorem 1.1 

We now prove Theorem 1.1, which will be an easy consequence of Strichartz esti- 
mates and the dispersive inequality. Let e > be a small constant to be chosen 
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later, and let R 3> 1 be a large number depending on e to be chosen later. From 
(22) we have 

11^(0)11^(^3) <1- 

Thus (by monotone convergence) we may find a time T > e~ 2 > (depending on 
it) such that 

11^(0)11^(^^)^3) <e. (32) 
By hypothesis, we may then find a time T\ > T such that 



\u{T u x)\ 2 dx<e. 

IB(0,R) 

Then of course we have 



1 



mHT^x)^ dx<e, (33) 



where rjji was defined in the notation section. We now bootstrap this control of the 
local mass at a fixed time T\ to control of the local mass at nearby times. From 
(1) we observe the mass flux identity 

<9 t M 2 = -2V-Im(VmZ) 

and hence by integration by parts 19 



dt J VrH 2 dx = 2 / Vr) R ■ Im(Vuu). 

Since ||V?7fl||L°o = 0(1/R), we thus obtain from (3) the crude estimate 

\dt f VRHt)\ 2 dx\ < 1/R 

for all times t (far better estimates are available, but this is already sufficient for 
our purposes). Thus by (33) we see that 



sup / rj R \u(t,x)\ 2 dx < e 
*e[T 1 - £ - 1 /4,Ti] J 



if R is large enough depending on e. In particular we see that for any t E [Ti — 
£-1/4^ j^j we nave 

||«(t)IU3( B( o,*/2)) < ll«WII^ 2 ( R 3)hWII^ (B( o^/ 2 )) <e 1/2 
(using (28)), while 

H-WIL3 (R 3 VB(0iii/2)) < ll^ll^^^^^ll-WII^) <R~ 2/3 

using (3), (5). Thus if R is large enough depending on e we have 

sup IMU R 3) <e 1/2 - (34) 

te[Tl _ e -l/4 iTl] L ^ 

We now claim the following smallness bound on the linear evolution (starting from 
Ti). 



^^One can justify this step by starting with smooth solutions u and then taking limits using 
the H 1 local well posedness theory. We will justify similar formal manipulations later in this 
paper without further comment. 
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Lemma 4.1. We have the estimate 

LtL«([T ll + oo)xR 3 ) 



(t - Tl)A -(^)ll r , 4r ,^,_wR3^- 1/32 - (35) 



Proof We first apply Duhamel's formula (27) to obtain 



u(Ti) = e mA u(0) -if' e l ^-^ A F{u{t')) dt' 
Jo 



and hence 



i / 1 e^-^ A F(u(t')) dt' 

JTi-£- 1 /4 

i(*-*')A F ( u (t/)) eft'. 



(36) 



o 



From (32) we see that the contribution of the first term (the linear component) of 
(36) is acceptable. Now consider the second term (the "recent past" contribution). 
By Strichartz (22) and Minkowski's inequality we see that the contribution of this 
term to (35) is bounded by 



< 



r mum K s dt>. 

JTi-e- 1 / 4 



But by the Leibnitz rule and Holder we can bound this by 

£ \\ U hrLl\\u\\ L 2 L o.\\u\\ L 2 w l,6 

where all norms are over [Ti - e- 1 / 4 , Ti] x R 3 . But by Lemma 3.2 we sec that 

ll«lli?L»,||V«|| L 2 i2 <e- 1/8 
and so by combining this with (34) we see that this contribution is acceptable. 

It remains to prove that 

" ' e'<'-''» 4 FK0)<!*'ll,.„ m _w n .,S« 1/32 . 



Jo 



l L 4 L6([T 1 , + oo)xR ) 

On one hand, observe from Duhamel's formula that 

/ 1 e^-^ A F(u(t')) dt' = e 4 (*- Tl+e " /4 ) A ( U (T 1 - e- 1 / 4 ) - u(0)) 

Jo 

and so from Strichartz (20), Holder, and (3) we have 

II / 1 ' e^ A F(u(t')) *>\\ Lt xJV) < h(T 1 - £ - 1 /4 ) _ u(0) || i2(R3) < L 

J 

Thus by Holder's inequality again, it will suffice to show that 

II / " ^("(0) ^IL 4 ^ ([ T 1;+ oo)xR 3 ) S £l/16 - 

J 



ASYMPTOTIC BEHAVIOR OF FOCUSING NLS 



19 



But from the dispersive inequality (19) we have 

II jf 1 e^-^KO^'ILoo^ <l 1 " (t-tr" /2 \\F(u(t'))\\ LHR ^dt'. 

But by (28) we sec that ||F(u(i'))|| Ll(R 3 = 0(1), and thus 

II l 1 ' e^ A F{u{t')) dt'lLoo^j < (e- 1/4 + (t- T!))- 1/2 . 
The claim follows. ■ 

We can pass from the bounds (35) on the linear solution to a corresponding bound 
for the non-linear solution by standard arguments. Indeed, from the Duhamcl 
formula (27) we have 

u(t) = e i (*- Tl ) A «(T 1 ) - i f e^-^ A F{u{t')) dt', 

and upon taking LjLldTi, +oo) x R 3 ) norms, we see from (35) and Strichartz (23) 
that 

ll U llLfL S ([T 1 ,+oo)xR S ) ~ ^ + \\ F M\\ltWl-*'*([T u +oo)xR a y 

But by Leibnitz and Holder we have 

ll f, ( u )lli?w. 1 - 8 / 5 ([T 1 ,+oo)xR 8 ) ~ ''"HifLSdTx.+oojxR^I^Hirffidri.+oojxR 3 )- 

Thus from (3) and a standard continuity argument (e.g. replacing [Ti,+oo) by 
[Ti, T) and then letting T — > +oo) we thus see that 

IML^ar^xR 3 )^ 1732 

if e was chosen sufficiently small. Finally, we set u + as 

u+ := e- 4TlA u(T 1 ) - i / e - u ' A F(u(t')) dt'; 



the Strichartz inequality (23) shows that the latter integral converges in H 1 since we 
have just shown F(u) to lie in L 2 t wi' 6 ^ 5 ([Ti, +oo) x R 3 ). Then from the Duhamel 
formula again we see that 



/+oo 
e- u ' A F{u{t')) dt' 

for all t > Ti, and thus (by the convergence of the integral) we have 

ihnjKt)- e ^ + || ffl(R3) = o 



1 1 1 1 1 i a \ 1 1 — ( 11 A 

as desired. 



Remark. This argument can mostly be extended to the non-radial case, and shows 
in fact that one has scattering to a free solution whenever 

lim inf sup / \u(t, x)\ 2 dx < e 

for some specific e, R > depending on the energy E; in other words, one must 
always have some concentration of the L 2 mass in time in order to prevent the 
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solution from scattering to a free solution, although in the non-radial case one does 
not have much control on the location xq = xo(t) of this concentration (though 
the sub-critical nature of this equation, combined with finite speed of propagation 
heuristics, suggest that we can make x (t) depend in a Lipschitz manner in time). 
We leave the details of this generalization to the reader (similar results also occur 
in e.g. [29], [6]). 



5. Proof of Theorem 1.2 



We now prove Theorem 1.2, which is another application of Strichartz and disper- 
sive inequalities. We begin with a standard preliminary estimate, which in some 
sense asserts the existence of an adjoint wave operator. 

Lemma 5.1. Let u be any solution to (1) obeying (3) (not necessarily spherically 
symmetric), and let f G if _1 (R ). Then the limit 

lirn (u(t),e HA f) L 2 (37) 

t — >+oo 

exists and is bounded in magnitude by 0{\\f\\ H _ 1( ^s^). 
Proof It is clear that 

sup |<u(t),e itA /M<||/|| R8 

t€[0,+oo) k ' 

just from (3) and duality of i? 1 (R 3 ) and _ff _1 (R 3 ). Thus by the linearity of 
(u(t), e ttA f) L 2 in / and a standard limiting argument, it will suffice to prove the 
convergence of the limit (37) when / is a test function, i.e. a smooth, compactly 
supported function. In particular / now lies in L 1 and so from (19) we have the 
bound 

||e*' A /|U=o < C f t- 3 / 2 (38) 
for all t > 0, where Cj is a constant depending on /. 

It will suffice to show that (37) is a Cauchy sequence. Accordingly, we pick an e > 
and seek to find a large time T > such that 

|<u(t),e <tA /> ia - (u(T),e* TA f) L2 \ <e (39) 

for all t >T. To show this we use Duhamel's formula (27) to write 

u(t) = e^' T ^u{T) - i J' e^-^ A F{u{t')) dt' 

and hence the left-hand side of (39) can be rewritten as 

\(j\e i(t -^ A F(u(t<)),e* tA f) L2 dt'\. 

Thus it will suffice to show that 

r + oo 

J t \(F(u(t')) 7 e^ A f) L2 \dt'<e. 
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But from (28) we see that \\F(u)\\ L1{ - R 3 ) < 1, while from (19) we have ||e it ' A /|| ioc(R3) < 

(0~ 3/2 U/ll L i(R 3 )- Thc claim thus follows from the decay of (i')~ 3/2 if T is taken 
sufficiently large. ■ 

The expression (37) is thus a bounded linear functional (of /) on i7 _1 (R 3 ), and 
hence by duality there exists a unique u + G i/ 1 (R 3 ) such that 

km (u(t),e itA f) L 2 = (u + J) L 2 
for all / £ H~ , and we have the bound II^+II^i^r 3 ) ^ 1- If we then define u w ^ 

by 

M wbW ■= u{t) - e ltA u+ 
then we clearly have (6) and (8) (and hence (9)), thanks to the unitarity of e ltA . 
Also we obtain (7) from (3). Note that this argument also shows that u + (and 
hence u wb ) are unique. Also it is clear from construction that if u is spherically 
symmetric then u + and it wb are also. 

Now we prove (10). Fix r. Since 

(id t + A)u = F(u) 

and 

{id t + Ay tA u+ = 

we have 

(id t + A)« wb - F(u wh ) = F(u) - F(u wh ). 
Thus it will suffice to show that 

r Iim \\F(u) - F(u wh )\\ LlH1J[TT+T]xR3) = 0. 

Let us just show this for thc homogeneous norm H\, the contribution of thc L 2 X 
norm being similar. We observe the pointwise estimate 

\W(F(u)-F(u wh ))\ < (|Vu wb | + |Vu|)(|«|+|u wb |)|«-u wb |+(|«|+|u wb |) 2 |V(^-«wb)l 
and hence by Holder 

\\F(u) - F(u wh )\\ LlAl <|||Vu wb | + |Vw||| LrL 2||| u | + \u wh \\\ L 2 Lr \\u-u wh \\ L 2 Lr 

+ IIM + l«wbllli?iSlH M l + \ u wb\W LrL°\\V(u - u wh )\\ L 2 L 6 
where all norms are in [T,T + t] x R 3 . From (21) we have 

l|V(M-U wb )|| L?LS(RxR 3 ) <1 

and hence 

T lim ||V(u - «wb)IL ?iS([ T,T + .]xR 3 ) = °' 
A similar argument using (22) gives 

r Iim ||U - "wbllL ? L»([T,T+r]xR 3 ) = °" 

Meanwhile from (21), (22) and Lemma 3.2 we have 

IIM + l U wblllL ? LP([T,T+r]xR 3 ) ^ ^ 
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for p = 6, oo. Combining all these bounds together with (7) and Sobolev embedding 
we thus obtain 

T lim \\F(u) - F(u wh )\\ Llt/llJlTT+T]xR3) = 

as desired. The analogous argument for the L 2 X norm is similar (because every 
bound that we have on derivatives of u and u wb , we also have on u and u wb 
themselves) and is omitted. This proves (10). 

We now prove (11). Observe that 
^K H R 3 ) = ll"wb(*)l£ 2( R«) + H^+H^R 3 ) + 2Rc( Uwb (t), e riA U+ ) i2(R3) . 

The left-hand side is MASS(u(0)) by conservation of mass. The second term on 
the right-hand side is MASS(u + ) since e ltA is unitary. The third term on the right 
goes to zero as t — > +oo by (8), and the claim follows. 

Now we prove (12), which is in the same spirit but a little trickier. We begin with 
the H 1 version of the above identity, 

H Vu WI^(R 3 ) = l|V Uwb W|!^ (R 3 ) + ||V e 4 *V||2 2(R3) +2Re(V Uwb (t),V e lt V} i2(R 3 ) . 

Observe that the third term on the right still goes to zero by (8) (moving both 
derivatives over to e ltA u + ). The left-hand side is equal to 2ENERGY(u(t)) + 
\ J \u(t, x)\ 4 dx, while the first term on the right is equal to 2ENERGY(u w ^(t)) + 
\ I l ?i wb^' x )\ 4 dx. Finally, the second term on the right-hand side is just J |Vw + (a;)| 2 dx. 
Putting all this together we see that (12) will follow if we can show 

lim / \u(t,x)\ 4 - \u wh (t,x)\ 4 dx = 0. 
t ^+ oc jR 3 

From (3), (7) we see that || U (*)|| L 4(R 3 )) ll u wbllL 4 (R 3 ) ~ ^' Thus the claim follows 
from the pointwise bound 

\\u(t,x)\ 4 -\u wh (t,x)\ 4 \<(\u(t,x)\ + \u wh (t,x)\) 3 (\u-u wh \(t,x)) 

(40) 

= (\u(t,x)\ + \u wh (t,x)\f\e* tA u + (x)\ (41) 
and Lemma 3.3 (applied with p — 4). ■ 



Remark. One can obtain more explicit Duhamel-style formulae for u wb and u + , 
namely 



/i-oo 
e i(t-t')A F{u{tl)) df , 



and 

/■+oo 

u+ = u(0) - i / e- itA F(u(t)) dt 
Jo 

although these integrals do not converge in the energy class, but only in weaker 
senses such as in the distributional sense (which is what is essentially done above). 
Note that the Duhamel-type formula for u wb consists entirely of a forcing term 
(from the future) and no linear term; we thus expect w wb to be smoother and 



p + OC 

t 

time decay of F(u) does cause some difficulty in making this heuristic precise. As 



more localized than u itself, although the infinite time integral J t °° and the lack of 
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indicated in Theorem 1.3, we will throw away a small energy error from u w ] J (which 
is basically caused by the coupling of the bound state with either the initial data 
e ltA u(0) or the data at infinity e ttA u + , as well as the non-linear self-interactions 
of the radiation term e ltA u + ) before establishing our final smoothness and decay 
estimates. 

Remark. The above result establishes L 2 and H 1 type estimates on the radiation 
term u+. It is an interesting problem as to whether there are similar estimates for 
higher regularities; in particular, if the map u(0) i— > u+ is bounded from £P(R 3 ) 
to H S {J\?) for some s > 1. This type of result is true in the defocusing case (see 
e.g. [15]), but our methods do not seem strong enough to establish this result for 
the focusing equation. The obstruction is a scenario in which the bound state (w w ] 3 
or Ufo) from continually emits minute amounts of high frequency radiation (which 
could conceivably cause the H s norm to become unbounded) without contradicting 
conservation of mass and energy, and without causing the bound state to collapse 
into the vacuum state. We do not know how to show that such a scenario does not 
occur. 



6. Decomposition into incoming and outgoing waves 

We now introduce a fundamental concept in our analysis of the free Schrodingcr 
propagator e ttA , namely the decomposition of a spherically symmetric function 
into incoming and outgoing components (plus a smooth error) ; these basically cor- 
respond to superpositions of radial waves e 2 "' 2 ' with u> < and u> > respectively. 
A similar decomposition (constructed using pseudo-differential operators) is also 
available for non-spherically-symmetric functions, see for instance [35]. 

Proposition 6.1. Let 5 > 0, and let R > be such that either R = or R > 1. 
Let f be a spherically symmetric test function supported on the exterior region 
{x : \x\ > R} (so if R = then f is just any spherically symmetric test function). 
Then there exists a decomposition 

f = /+ + /- + f smooth (42) 

with the following properties: 

• (L 2 -boundedness) We have 

ll/+ll L2( R»). H/-lli2(R 8 ). WfsmoothW^^ < H/lli2 ( R8)- ( 43 ) 

• (Smooth error) The function f smooth obeys the estimates 

||V/ smoot ,|| L2(R 3 ) < {Rf-'^WfW^^ (44) 
for all a > 0. 

• (Asymptotic decay of incoming waves) For any u e L 2 (R 3 ) and any e > 0, 
there exists a time T = T(u , e, S, R) > such that 

\(e itA u ,f.)\<e\\f\\ L2(R3) (45) 
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for all t > T and all test functions f on {x : \x\ > R}. (Informally, this 
means that as t — > +00, the free wave e 4 * A uo consists asymptotically of 
purely outgoing radiation). 
• (Exponential decay near origin in the favorable time direction) We have 

° \W^ A f + \\ L2{B{0 , R/8)) dt < (R)-"\\f\\ H ^ (W) (46) 



1 1 

and 

f° WVe^-^f-WmB^R/,)) dt<{R)-*\\f\\ ^ (47) 

J — oo 

for any a, (3 > 0. (Thus, outgoing waves radiate away from the origin as 
t — > +oo, while incoming waves radiate away from the origin as t — > — oo J 
• (Local smoothing in the favorable time direction) We have 

/+oo 
\\^r 2 ^ to)A u\\ H -^) dt s (R)- 1+5 \\f\\ H -. r+S{W) (48) 

and similarly 

f_l ll^)- 2 e <(t - t0)A /-ll^ )( R 3 ) * S W-^ll/H^-^cR-) ( 49 ) 



/or any a > 0. ffere i?#(R ) is the scaled Sobolev space defined in the 
notation section. 



In the above estimates we allow implicit constants to depend on the exponents S, a, 
(3 but not on R. 



Proof In what follows the reader may find it helpful to keep in mind the following 
heuristics: firstly, that the Littlewood-Paley projection operators Pj only introduce 
a spatial uncertainty of 0(2~ J ), by the uncertainty principle; and secondly, the 
dispersion relation indicates that the Schrodinger propagator e 4 (*-*o)A^ wnen acting 
on functions of frequency ~ 2 3 (such as Pjf), will move this function by speeds 
comparable to 2 3 . Because all functions here will be spherically symmetric, there 
are only two directions for propagation: inwards toward the origin, and outwards 
away from the origin. To separate these two modes 20 we will use (truncated) radial 
Riesz projections. 

In the argument that follows we assume some familiarity with the heuristics of 
stationary phase, and in particular the ability to get arbitrary decay 21 of oscillatory 

2t ^It is possible instead to proceed using the full machinery of Fourier integral operators and 
semi-classical analysis; for instance, f± is basically the projection of / to the phase space region 
51 ^ {x)~ 1+s , ±5 ■ x > 0, while f smoot h i s basically the projection to the region |£| < {x)~ 1+s . 
However, to keep the paper reasonably self-contained, we have proceeded more explicitly (though 
perhaps at the cost of brevity), relying mostly on integration by parts, stationary phase heuristics, 
and the Fourier inversion formula. The reader who is content to accept this Proposition may in 
fact skip the technicalities and move on to the proof of Theorem 1.3 in the next section. 

21 We shall refer to this as the principle of non- stationary phase. We also need a Van der 
Corput lemma, which asserts that Jpj3 e 11 ^^ ip(x) = 0(A — 3 / 2 ) if <j> has exactly one stationary 
point with V 2 > A at that stationary point, i/j is a rescaled bump function, and <p obeys the 
usual smoothness bounds on the domain of ip. See [42]. 
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integrals J e l ^ x 'ip(x) dx via repeated integration by parts when the oscillation 
of the phase exceeds the amount one loses when differentiating ip; we refer to [42] 
for more precise formulations of this heuristic. 

Let us first consider the case when R = 0, and / is supported on B(0, 2). In this 
case we set /_ := f smooth ■= and /+ := /. The bounds (43), (45), (44), (49) 
are then trivial, while (46), (47) are vacuously true (since P(0, P/8) is empty). It 
thus remains only to verify the outgoing local smoothing estimate (48). We have 
to show that 

IK*>-V<«-*>*/|| R , dt<\\f\\ H _ a _ 1+5{K s y 

I to 

We use Littlewood-Paley operators to split / = P<of + J2j>o ^j/' where P<of ■= 
J2j<oPjf- Consider first the contribution of the low frequencies P<of- Observe 
from (19) that 

||e i( *- to)A P<o/||oo < \t t \-^\\P< f\\ LHB{0A)) 

<\t-t \- 3 / 2 \\P< a f\\ L 2 

<\t-to\- 3/2 \\f\\ H -«-^, 

and also by Bernstein's inequality 

||e i(t - to)A P<o/||oo < ||e 4( *- to)A P<o/||L 2 = \\P<ofh' < \\fWH-~-w 
and hence 

||e*(*-*o)Ap< o/ || oo < (t-t Q )- 3 / 2 \\f\\ H - a -i +s . 
By Holder we thus see that 

\\(x)- 2 e^-^ A P< f\\ L2 < (t-t )- 3 / 2 \\f\\ H - a - 1+ s 

and the claim follows by integrating in t (noting that L 2 controls H~ a ). To control 
the high frequencies, it will suffice to show that 

r+oo 

/ ||(x)- 2 e ^- t ") A P,/|| H _ Q(R3) dt < 2(-«+WW||p./|| i2(R3) < 2- 5 -'/ 2 ||/|| H _ a+1 _, (R 3 ) 

J to 

for all j > 0, since the claim then follows by summing in j. Fix j; we define 
t = Tj := 2^ 1+s / 2 ^ and divide the above integral into the immediate future J t *° +T 
and the later future r + °° . For the immediate future, we use Lemma 2.1 to discard 
(x)~ 2 , and thus bound this contribution by 

rto+T 

dt, 



< / " \\e lit - to)A Pif\\ H ^ 
Jto y 



which is then acceptable by definition of r since e'(* io ' A preserves the H a norm. 
Now consider the later future J^+ T - For this contribution we write Pj = PjPj where 
Pj is a slight enlargment of the Fourier projection operator Pj. Now observe that 
using stationary phase that the convolution kernel K(x) of e *( t -*o)A)p^ j s mos ^iy 
concentrated in the region \x | > 2i\t-t \ > 2«/ 2 (the latter bound following 
definition of t); indeed from the principle of non-stationary phase we can obtain 
bounds of the form \K(x)\ < (2 J |i - t \)- M in the region \x\ < 2' J \t - t \. Thus the 
contribution of the region of space |x| -C 2 J \t — to] is easily seen to be manageable. 
Thus by applying a smooth cutoff we may restrict to the region \x\ > 2?\t — to | ■ 
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But then we may invoke Lemma 2.1 again to replace (x) 2 by (2 J \t — to\) 2 ■ Thus 
we can bound this contribution by 

+T (2 3 |t-to|)- 2 || e 4( *- to)A P,.f|| ff _„ (R 3 ) ^ 

which by the boundcdness of e l (*-*o) A j s bounded by 

2- 2j r- 1 ||P,/|| H _ a(R 3 ) 

which is acceptable by definition of r. This concludes the treatment when / is 
supported on B(Q, 2). 

It thus remains to consider the case R > 1, since the case R = then follows from 
the -R > 1 case and the localized case just considered. We can then replace (R) by 
the comparable quantity R throughout. We then define f smooth as 

} 'smooth ■= ^ P 3-f + r lR E P if> 

where t]r is defined in the notation section, and now proceed to verify the bounds 
(44) on f smooth- The contribution of ^2 2 j<R s - 1 Pjf 1S cl ear just from taking Fourier 
transforms, so it will suffice to verify that 

E HWi/ii^R.) ;$tf (5 - 1)(1+a) ll/ll*-»(Rv 

But / is supported on R 3 \_B(0, R), which is a distance R/2 from the support of 
7]r, while the convolution kernel of Pj decays rapidly outside of B(0, 2 - - 7 ), as does 
all of its derivatives, so it is easy to verify the pointwise estimate 

IIV^/IU < (VR)- M \\f\\ H _ a{R3) 

for any M > 0. Summing over all j with 2- 7 > we obtain the result (if M 

is sufficiently large depending on 8, a we may absorb all powers of R that arise on 
the right-hand side). This proves (44). 

It then remains to decompose 

(i- VR ) £ Pjf = f- + U 

2i>R f >- 1 

with the desired properties. We shall do this by decomposing each (1 — f]R)Pjf 
separately for each 2 J > R 5 ~ x , 

(1 - m )P 3 f = /_,,■ + /+,, (50) 
and then defining f± := Y,2i>RS-i f±,j- 

We begin by taking the radial Fourier transform gj of the odd extension of (1 — 
VR)Pjf, 

9j (p) := / + °° e - 2 ^sgn(r)(l - m(\r\))Pj f(\r\) dr, (51) 

J — OO 
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so that gj is an odd function on R, and observing the Fourier inversion formula 

/oo 
gj(p)e p dp 
-oo 

on R 3 , where e p is the radial function 

e p {x) := e 2nip W. 

We can thus define 

f±,j ■= 9j(p){l -VR/2)e p dp (52) 

J±[0,+oo) 

and we thus obtain the decomposition (50) (and hence (42)); note that the factor 
(1 — VR/2) does not affect (50) since it equals I on the support of (1 — r]n). 

Now we verify the L 2 bounds (43). It suffices to show the claim for f±. By (52), it 
suffices to show that 

II j + ^ - VR/2K d / 3 ll L 2(R 3 ) < ll/ll L 2(R3) (53) 

where g :— ^2i2j>r s - 1 9j- We discard the (1 — f]Ri<i) and use polar co-ordinates to 
write the left-hand side as 

< ||r / g(p)e 2mrp dp\\ L , 
Since g is odd, we can integrate by parts to estimate this as 

< II / 9'{p)e 2mrp dp\\ L 2 

J±\0,oo) 



'±[0,00) 

which by Plancherel is bounded by 

<(J K W(P)\ 2 dp) 1 ' 2 . 
Expanding out g' using (51), this is 

< ( / I / e - 2 -"sgn(r)(I - m (\r\)) E P i/(M) rdr f 

JR. J -00 2J>_R' S - 1 

which by Plancherel again is bounded by 

/+00 
1(1 ~VR(\r\)) E Pif{\r\)\ 2 r 2 drY' 2 . 
-00 2i>R s - 1 

But if we discard the cutoff (1 — rm) and undo the polar co-ordinates, this can be 
estimated as 

11 E pjfWw) z 11/11^ 

2i>R 8 - 1 

as desired. 

Now we prove the asymptotic decay estimate (45) for the incoming wave. From 
(43) we already have 

supKe^uo,/-)^^)! < ll/ll L 2 ( R3 ) ||wo|| L2( R^ ) , 
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so by the usual Ricmann-Lcbesgue type limiting argument (as in Lemma 3.3) it will 
suffice to prove (45) assuming that uo is a (spherically symmetric) test function. 
But in that case e ltA uo has the asymptotics 

e ltA u (x) = t-WjW^aix/t) + o L2(R 3 } (l) 

where a is a spherically symmetric Schwartz function (it is essentially the Fourier 
transform of uo) and ° i 2(R j 3 )(l) i s an error whose L 2 (R 3 ) norm goes to 1 as t — > 
+cxo; sec e.g. [43] (or one can work directly from (18)). Using (43) again, it thus 
suffices to show that 

t -3/2| (e i|x|V4* o(a . /t))/ _ )i2(R8) | < £ ||/|| L2(R3) 

for all / e L 2 (R 3 ) supported on R 3 \5(0,i?), and all t > T = T(a,s,R). We 
expand out /_ using (52) and estimate the left-hand side by 

<t- 3/2 f \9{p)W^ 2/U a{x/t),{l-r lR/2 )e p ) w \ dp. 

J — OO 

We can write this in polar co-ordinates as 

<t- 3/2 f W)\\ I e^/^ r 2 a[r/t){l _ 1lR/2{r)) dr \ dpi 
and then break this up dyadically as 

oo -o r 

Z £*~ 3/2 / \9(P)W / e" 2 / 4t -^r 2 a(rA)( %fcK - %*-i fl (r)) dr\ dp, 
k=o J -°° J (54) 

But the phase ir 2 /At — ipr oscillates in r with frequency at least \p\ + 2 k R/t on 
the support of rj 2 k R — f] 2 k-\ R (because p is negative), and so by the principle of 
non-stationary phase we have the bounds 



,ir 2 /4t-ipr 2 



r 2 a(r/t)( V2 , R ^ V2k - lR (r)) dr\ < (2 k \p\R + 2 2k R 2 /ty M (2 k R) 3 



for any M > by repeated integration by parts (the function a(r/t) obeys symbol 
estimates in r regardless of what t is, because a is Schwartz). This decays rapidly 
for \p\ ^> R~ 1 2~ k . Hence we can estimate (54) as 

oo 

<^i- 3 / 2 (2 2fc i? 2 A)" M (2 fe i?) 3 (^ 1 2-' £ ) 1/2 ||.g|U ? . 

This decays rapidly for 2 k > i 1 ' 2 ./? -1 , so we may sum in k to obtain the bound 

< r^it^R^RfiR-h-^R^WgU* = t- 1/4 ||.g|| L 2. 
But by (51) and (one-dimensional) Plancherel we have the somewhat crude estimate 
\\9ihi < 11(1 - m)P-jf\\ui < \\Pjf\\ L3(R?) (55) 

and more generally 

NU ? <||(i-rm) £ p jfh^W E p jf WmW) Z H/H^r 3 ) 

and the claim follows for t sufficiently large. 
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Now we prove (46), which is basically an application of the principle of non- 
stationary phase. It suffices to show that 

r- + oo 

|[Ve ,( t -to)A /+j||i2(B(oi?/8)) dt<{VR)- M \\P 3 .f\\ L2(R3) 

to 

for every Z> > R 6 - 1 and every M > 0, since the claim then follows by taking M 
large enough depending on a, (3,5 and summing in j. 

Fix j. By time translation invariance we may take to = 0. We expand the left-hand 
side using (52) and use Minkowski's inequality to estimate by 

r- + co r- + co 

< / lft(p)ll|Ve i * A (l-? ? i J /2)ep|U2 (s(0iJ j /8)) dpdt. 

Jo Jo 

It will suffice to prove the estimate when (1 — 7]r/2) is replaced by (rju — 77^/2), since 
the claim will then follow by replacing R by 2 k R for k = 0, 1, . . . and summing the 
telescoping series. By (18), we have 

e ltA (m VR/2)e P (x) = ^- J e^-y\ 2 '^{ m - V R /2(y))e 2vipM dy; 

using polar co-ordinates in y this becomes 

W/** / e-^ rcos6 / 2t ( m - VR/2 (r))e lr2 / 4t e 2m P r sin9d9r 2 dr. 
Jo Jo 



W 2 ^ 



Substituting s = cos 9, this becomes 

S I f-1 poo 

° J\x\ 2 /4t / / ir 2 /4t p -i\x\rs/2t 1-Kipr, 



t 3/2 



/I f-OO 
J e ^ 2 /4t e ^\x\rs/2t e 2 m pr^ R _ 7?fl/2 ( r )) ^ d rds. 



If x £ B(0, R/8), then the phase e ir2 /*t e -i\x\rs/2t e 2*i P r osci u ates m r at a rate of at 
least > p+R/t on the support of 777?— f? fi / 2 (note that this uses the positivity of p, i.e. 
the fact that we are only considering outgoing waves). Meanwhile, every derivative 
applied to (rjn — rj R / 2 {r))r 2 lowers this quantity by a factor of about R. Thus one 
can integrate by parts repeatedly to obtain a bound of 0(t~ 3 ^ 2 R 3 (pR + R 2 /t)~ M ) 
for any M'. Applying \7 X introduces factors of the order of at worst R/t and thus 
we obtain a bound of the form 

\Ve ltA ( m - V R/2)e P (x)\ < t~ 5 / 2 R 4 (pR + R 2 /ty M . 

Thus it remains to show that 

z' + oo r+oc 

J J t-^ 2 R\ P R + R 2 /t)- M '\g {p)\ dpdt < (2^)- M ||P J /|| i2(R3) . 
We can evaluate the t integral to estimate the left-hand side as 

r+oc r>2 

l { W) r3/2Ri{pR) ~ M lgM dp ' 

If p was much larger than 1/R, e.g. p > (2^/i?) 1 / 2 , then the claim is now easily 
proven from (55), since 2^ > R s ~ x . So we may restrict our attention to the region 
p < (2 J /i?) 1 / 2 . But then we may use (51) to write 

/ + OC 
e - 2 ^sgn(r)(l-r ?fl (|r|))P ) /(|r|) dr, 
-OO 
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which by collapsing the integral to J °° and then undoing the polar co-ordinates 
becomes 

9iiP) = C | r3 s^\x\p) l -^ 2 {x) P 3 f(x) dx. 

We now integrate by parts repeatedly, integrating Pjf to effectively gain a factor 
of 2~- 7 (by the frequency localization) while differentiating sin(27r|x|/9) or ^—^r^-, 

at worst losing a factor of (2? /R) 1 / 2 . This will eventually gain us an arbitrary 
power of (2 : >R)~ M , and then if we estimate everything else crudely (e.g. using 
Cauchy-Schwarz and Plancherel to handle the p integral) we will eventually show 
this contribution is acceptable (recall that 2 3 i? > R s and so we can eventually 
overcome any polynomial losses in R). This proves (46). The proof of (47) is very 
similar and is left to the reader (it is essentially the conjugate of (46)), so we now 
turn to (48) and (49). Actually it will suffice to prove (48) since (49) is essen- 
tially the conjugate of (48) (recall that conjugation turns the forward Schrodingcr 
evolution into the backward Schrodingcr evolution). 

We now prove (48). Again we may use time translation invariance to set to = 0. 
First observe that the portion in the ball B(0, R/8) can be easily treated by the 
estimate (46) just proven, so we may freely insert a cutoff (1 — i] R / s ). By dyadic 
decomposition in j it will suffice to show that 

r+oo 

J a 11(1 - VR/s)(x)- 2 e itA f + J H . a{R3) dt < J?- 1 + 6 ||P j /|| /f _ a _ 1(R 3 ) 
for all 23 >R s -\ 

Fix j. The intuition here is that at time t, the wave f+j will have propagated away 
from the origin by a distance of ~ 2H, and will hence the weight (x) -2 will have 
size 0((2H + R)~ 2 ). Integrating this in time gives us Rr 1 2~i , which explains both 
the gain of the power of R and the smoothing of one derivative. This should be 
contrasted with the standard Kato local smoothing result [38], [48], which works 
both forward and backwards in time (for both incoming and outgoing waves) and 
gives square-integrable bounds in time, but gains only half a derivative and no gain 
in R. 

We turn to the details. We first split the time integral / + °° into the near future 
J T and the distant future f^°, where 

r := 2-iR 1+s ' 2 (56) 

(this can be interpreted as the time it takes for the wave to leave the region |x| <~ R, 
times a safety margin of i? 5 / 2 ). We first consider the near future. For this region wc 
estimate (1 — Vr/8){ x )~ 2 by R~ 2 (using Lemma 2.1) and estimate this contribution 

J\- 2 \\e^f + J HRa{K3) dt. 
Since e ltA preserves the H^" norm, it will thus suffice by (56) to show that 

ll/^H^R 3 ) £ ll^/II^R 3 ) ( 57 ) 
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(this estimate will also be useful for the treatment of the distant future, below). 
We expand this using (52) as 



<?»(! -VR/2)e p dp\\ H - a(B s < 2-^1^/11 R s 

[0,+oc) { ' 

By an easy modification of the proof of (43) we already have that 

~VR/2)e p dp\\ L2R s < ll^/ll L 2 f R» v 
'[0,+oc) y ' y 1 

so by Littlewood-Paley decomposition (and the hypothesis V > i? 1_(5 ) it will suffice 
to prove the estimates 

II E P r I 9j{p)0- - VR/2)e P dp\\ R3 < (2^)- M ||P,-/ 1| R3 
j'<j-io J [°' +oc ) (58) 

for all M; here of course the implicit constant is allowed to depend on M. 

This estimate is based on two observations; firstly, that the Fourier transform ((1 — 
VR/2) e p) is concentrated near the sphere |£| ~ p (and in particular Py ((1 — rj R / 2 ) e P ) 
is very small if p ^> 2 J ), and secondly that gj{p) decays very quickly if p/2 3 is small. 
We begin by quantifying the first observation. We write the Fourier transform 
(1 — r lR/2) e p) h {€) °ut as a telescoping sum 

E L(^R(y)-V2^R(y))e 2m{plyl - i - y) dy. 

k=0 JIi 

If p > |£| + then the phase oscillates by at least p — |£|, and so by the principle 
of non-stationary phase we can bound this expression by 

< J2 0{{2 k Rf( p -^)- M ) = 0(R*(^\)- M ) 

k=0 

for any large M. In particular we see from Plancherel that 

II E ^'((l-^hOIL^ <{(V+p)R)- M 

j'<j -10 

when p > 2 3 ! ~ 5 for any M (note that this decay can absorb any bounded powers 
of 2 3 or R). This easily allows us to control the portion of (58) arising from the 
region p > 2- 7 ~ 5 , using (55). 

Now we consider the contribution where p < 2 ] ~ 5 . For this contribution we discard 
the projector X^/<j-io an d instead establish decay of gj(p). By (51), (17) and 
the Fourier inversion formula we have 

f r+oo 

9iiP) = / 3 fe +1 (0 - M0)/(0 / e- 2 ^sgn(r)(l - VR (\r\))e 2ni ^< drd£ 

where ei = (1,0,0) (say). We can divide the r integral into positive and negative 
axes, and then decompose further dyadically (note there is no singularity at the 
origin because of the cutoff 1 — i] R (\r\). If p < 2- 7 ~ 5 , then the phase oscillates at 
a rate of at least > 2 3 , because of the support of the Littlewood-Paley multiplier 
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7723+ 1 — V2j ■ Thus by integrating by parts as before, we may bound this expression 

by 

\gj(p)\ < / R3 - v»(o\\m)\R(VRr M 

for any M > 0. Using Holder and Plancherel, and using the (2 J i?)~ M factor to 
absorb any powers of R and V which appear, we thus obtain the bound 

\ 9j (p)\ < (2^)-^||P J /ll L2(R3) (59) 

for all p < 2 J . In fact a similar argument also gives similar bounds for arbitrarily 
many derivatives of gj (p) in this region. By another application of Plancherel (in 
the radial direction) we thus have 

f 2 j ~ 5 

»(P)(1 - VR/2)e P dp\\ L2(R3) < (2^)- M ||P j /|| i2(R s ) 

(again absorbing any powers of 2 J or R which appear into the (2 J i?) _M factor, and 
so this portion of the p integral is also acceptable. This proves (58), and concludes 
the treatment of the near future. 

Now we consider the distant future J T °° ; we have to show that 

/+00 
ll(i - VR/s)(x)- 2 ^f + J HRa(B?) dt < R-^mfW^-i^ 

We divide the time integral further into dyadic blocks; it will suffice to show that 

2 k + 1 T 

I 11(1 - VR/s)(x)- 2 e itA f + J H - R s dt < 2- Sk R- 1 + s \\P j f\\ H . a . 1R3 

J2 k T R R 

for each k > 0. 

Fix k. Recalling the heuristic that the wave f + j propagates outward at speed 2 J , 
we expect the function f + j to be concentrated in the region |£| > 2' J+k r 3> 2 k R. 
Accordingly, we split the cutoff (1 — r/ R / 8 ) into the near component rj 2 k R — r)R/g 
and the far component 1 — rj 2 k R . 

Let us deal with the far component first. We can estimate (1 — (x)~ 2 by 
(2 k R)~ 2 (using Lemma 2.1), and so we can control this contribution by 

< [ {2 k R)- 2 \\e^f + ,\\ H - R3 dt. 

J2 k T R 

Using the boundedness of e ltA on H^ a and (57) (and then (56)) we can bound this 

by 

< 2 k T(2 k R)-*2-^\\P j f\\ L2{R3) ~ 2- k R-^\\P J f\\ HR ^ 1{Il3) 
which is acceptable. 

Now consider the near component, which will turn out to be very small. We expand 
out f + j using (52) to estimate this contribution (using Minkowski's inequality) by 

2 k + 1 T 

< [ I MP)\\\{^R - VR/s){x)- 2 e UA ((l - VR/2)e P )\\ H - a(R?) dt. 

J2"t J[0,+oo) r > 
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We crudely estimate the H R a norm by the L 2 (R 3 ) norm, and use Cauchy Schwarz, 
combined with the observation that {r} 2 k R — Vr/s){x)~ 2 has bounded L 2 norm on 
B(0, 2 k R), to estimate the previous by 

<[ I \9Ap)\\W tA (^-VR/2)e P )\\ L ~ ( B{ a .2*R))dt; (60) 

we shall be able to get away with these very inefficient estimates because we will 
shortly gain a large power of (2 J i?)~ M which can absorb these losses. Let us 
first deal with those p for which p > 2 J_5 ; as we saw in the treatment of the 
immediate future, this is the dominant case. We now compute e ltA ((l — rj R / 2 )& P ){x) 
for x e L 2 {B(0,2 k R)), using (18), as 

e ltA ((l - VR/2)e p )(x) = Ct-*' 2 f (1- VR/2 (x))e 2 ^y\ +i \ x -y\ 2 / 4t dy; 

we may decompose this dyadically as 

Ct ~ V2 E / 3 - m™-iR(x))e 2 * ipM+ilx - yl2/4t dy. (61) 

Observe that the gradient of the phase 2irp\y\ + \x — y| 2 /4i in y is 

V y (2np\y\ + \x- y\ 2 /U) = (p + M) ». _ £. ( 62 ) 

Since p > 2^ 5 , a; G B(0,2 k R), \y\ ~ 2 m i?, and f - 2 k r > 2 fc i?/2 J ', we thus see 
that (62) has magnitude > p + 2 m R; the point is that x/2t has magnitude at most 
0{R/t) <C O(p) and thus does not significantly affect the magnitude of the (62). 
Thus by repeated integration by parts we may estimate 



|e JtA ((l - VR/2)e P )(x)\ < t-*' 2 E (2 m Rf(2 m RpY 



-M 
rn=0 

for any M; in particular we have 

|e JtA ((l^^ /2 ) ep )( 2 ;)|<(p J R)- M 2- 3 ' £ / 2 

(using the the fact that p > 2 J_5 and 2 J > R 8 ^ 1 to absorb all other factors using 
the (pR)~ M term). Combining this with (55) we see that this contribution to (60) 
is 

< 2 - fc/2jR --|| J p. / || jf __ 1(R3) 

for any M > 0, as desired (note that the factor of (pR)~ M can absorb any losses 
of 2 J or R which appear). 

Finally we consider the portion of the low frequencies where p < 2 J_5 . Here we will 
enjoy unlimited decay in the 2 3 and R parameters (because of (59)) so the only 
issue is to obtain the required decay in the 2 k parameter, which basically amounts 
to ensuring that one obtains the expected i -3 / 2 decay in the L°° norm. Again we 
expand e rtA ((l — Tjfi/2)e p ){x) as (61), but divide now into two cases, depending on 
whether m < fc+10 or m > fc+10. First suppose m > fc + 10; then the gradient (62) 
has magnitude at least > \y\/t ~ 2 m ~ k 2 j R~ s / 2 (using (56)), because the x/2t term 
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is much smaller than the other terms. Thus we may integrate by parts repeatedly 
and estimate this portion of (61) in magnitude by 

<t- 3 / 2 (2 m R) 3 (2 m R2 m - k VR- s / 2 )- M 

m>k+W 

for arbitrary M. Discarding the 2 m ~ k factor we thus see this expression is 0(t~ 3 l 2 2~ Mk (2 : > R)~ M ) 
for any M; inserting this bound back into (60) and using (59) we see that this con- 
tribution is acceptable. 

Now suppose that m < k + 10. Here the gradient (62) might vanish, but it can only 
do so at one point (when y is parallel to x and p + \y\/2t = \x\/2t). Furthermore 
we have the double derivative estimate 

X7l(2irp\y\ + \x-y\ 2 /4t)>l/2t 

where we use A > B to denote the statement that A — B is positive definite; this is 
basically due to the fact that the function 2-xp\y\ is (non-strictly) convex, because 
p is positive. Hence by Van der Corput's lemma (see e.g. [42]) we obtain a bound 
of 

I / ^R-V2^R(x))e 2 ^ M e^~y\ 2 ^ dy\<r 3 / 2 ; 

summing this over m < k + 10 we thus see that this contribution to (61) is at most 
0((k)t~ 3 / 2 ). (Actually we may easily get rid of the (fc), but it will not be important 
in our argument in any event). If we combine this with (59) we obtain the desired 
result. 

This proves (48). As mentioned earlier, the proof of (49) is similar, so the proof of 
Proposition 6.1 is now complete. ■ 



7. Exterior energy decay of «wb 



We now begin the proof of Theorem 1.3, which will be done in a number of stages. 
In this section we will establish the following preliminary decay property of u ^ 
away from the origin: 

Proposition 7.1. Let the notation and assumptions be as in Theorem 1.3, and let 
< S < 1 be arbitrary. Then for every R > we have 

km sup f \Vu wh (t,x)\ 2 dx<R- 2+2S . 

Note that this Proposition only gives decay of the energy of u w ^ away from the 
origin; we suspect that the mass of u w ^ also decays away from the origin (which 
would allow us to upgrade the o^i(l) error in (13) to the more natural o H i(l)), 
but we do not know how to do this. This Proposition is significantly weaker than 
Theorem 1.3 but is also somewhat easier to prove, and the basic strategy used in 
this argument will be used again when we prove the rest of Theorem 1.3. 
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Proof Fix 6, R. The claim is trivial from (7) if R < 1, so we may take R ^> 1. Let 
To > be a large time depending on R, S to be chosen later. We need to show that 

\Vu wh (t ,x)f dx<R- 2 + 2S 

R \B(0,R) 

for all to > To. By duality it suffices to show that 

\(Vu wh (to)J) L2(n3) \<R- 1+5 (63) 

for all to > To, and all test functions / supported outside of -6(0, R) with bounded 
norm m L 2 (R 3 ). Since u w fo is spherically symmetric we may take / to be also. 

Fix /, t - We use Proposition 6.1 to decompose / = /_ + /+ + f smooth, and deal 
with the contribution of each term separately. The term f smooth is easy, in fact we 
just use integration by parts and (7) to bound 

KVw wb (to),/ S moot/ t ) i 2(R 3 )l ^ l|V/ smoo th|| L2( R_3 ) 

which is acceptable by (44). 

Now we consider the contribution of the outgoing wave /+; here the strategy will 
be to compare u with the asymptotic solution e ltA u + . By (6) it suffices to show 
that 

\(Vu(t Q ),f + ) L2{R3) <V e ^ A u + ,./+) L2(R3) | < R-i+ s . 
To prove this we first apply (27) to obtain 

rT 

u(to) = e t(t "~ T ^ A u(T) + i / e t(t "~^ A F(u(t)) dt 

J t 

for any T > 0. Taking inner products with / + , we thus obtain 
(Vu(io), f+h* = (Vu(T),e^ T -^ A f + ) L , + i [ f (F(u(t)),e^-^ A f + ) L2 dt. 

Jto J 

As T — > +oo, we have 

lim (Wu(T),e^ T -^ A f + ) L , = {Vu + ,e~^ A f + ) L , = (We^ A u + , f + ) L , 

thanks to (6), (8). Thus it will suffice to show that 

(VF(u(t)),e^-^ A f + ) L2 \ dt<R- 1+s . 



to 

We now observe that 

on the ball B(0, 1) this just comes from (3), (28) and Leibnitz, while outside of this 
ball this comes from (5) and (3). The claim then follows from Cauchy-Schwarz and 
(48) (with a = 0). 

Finally, we consider the incoming wave /_ ; which is similar except that we compare 
u to the initial data e ttA uo instead of the asymptotic data. From (45) we have 



and 



|(v e ^ A u+ ,/-) i2(R3) l<^ 1+5 

\{v^u{v),n L2(R?) \<R- l+5 
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for all to > To, if To is chosen sufficiently large. Thus by (6) it will suffice to show 
that 

\(Vu(to),f-) La(Ra) (V e rf » A U (0),/_) L2(R3) | < R- 1+s . 
Applying (27) we thus reduce to showing that 



£ | J (\7F(u(t)),e^-^ A f„) L2 \ dt < R- 



i+s 



But this follows from (49) (with a = 0) by arguing as in the / + case. This concludes 
the proof of (63) and hence of the Proposition. ■ 



8. Improving the smoothness of m wb 



We now use an iteration argument to bootstrap the asymptotic regularity of u w ^ 
from H 1 to H 1+a for arbitrary values of a, and also to gain some symbol-type decay 
estimates on u w ^ away from the origin. For inductive purposes it is convenient to 
phrase this regularity in a dual formulation. 

Theorem 8.1. Let the notation and assumptions be as in Theorem 1.3, and let 
S > and a > be arbitrary. Then for every e > there exists a time T = 
T(u, 5, a, e) > such that one has 

\(Vu wh (t)J) L ,\ < (R)-< 1+a X 1 -»\\f\\ H -.+e\\f\\» (64) 

for all t > T, all R > and all test functions f supported on R 3 \£?(0, R). Here 
H(R)~ a * s ^ e scaled Sobolev space defined in the notation section. 

At the end of this section we will use this Theorem to construct and thus prove 
Theorem 1.3. 

Proof We fix 5 > and induct on a (keeping e, R free to vary). It suffices by 
interpolation to prove the claim when a is a non-negative multiple of 5. When a = 
the claim follows immediately from (7) (just using the first term on the right-hand 
side of (64). Now suppose inductively that a > 5, and the claim (64) has already 
been proven for a replaced by a — S, and all e, R (assuming T sufficiently large 
of course). The idea is to use local smoothing and dispersive type estimates and 
Duhamel's identity (27) to boost the regularity of u w ^ by at least 5 of a derivative 
(the point being that u w ^ does not contain the linear part of the solution e ltA u + 
as t — > +oo, which of course has no gain in regularity except in the local sense of 
Lemma 3.3). 

Now fix s, and let T > be a large time (depending on all the above parameters) 
to be chosen later. We first obtain some bounds on the non-linearity F(u) away 
from the origin. 

Lemma 8.2. Let R > 1. Then there exists a T > such that 

|(V((1 -VR/s)F(u(t))),f) L2 \ < R-( 1+a -W-V\\(x)- 2 f\\ a+s +s wo \\fhi 

(65) 
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Proof We first verify this in the case where / vanishes on the ball -8(0, £ 200 ). 
Then on the support of / we have u = O(e wo ) by (5), and hence by (3) we have 

||V((1 - r 7fl/8 )^(n(i)))|| i2(R 3 ) < £ 100 

and the claim follows. Hence we may assume without loss of generality that / is 
supported in 5(0, e~ 200 ). Indeed by dyadic decomposition it then suffices to show 
that 

I <V((^/4 - vn/*)F(u(t))), DlA < R-*R-V+°-W-6) \\f\\ +s + s 200 \\f\\ L , 

(66) 

for all 1 < R < e~ 200 equal to a power of 2, and all test functions /, as the claim 
then follows by summing using the triangle inequality (and taking T to be the 
supremum of all the individual values of T obtained from each of the finite values 
of R). 

Fix R. By the inductive hypothesis, we can ensure that 

| {Vu wh (t), f) L2 1 < jr \\f\\ H - a+ . + e 500 \\f\\ L 2 

for all t > T and all test functions / on the annulus B(0, R)\B(0, -R/64), if T is 
chosen large enough. By duality this implies that 

(VR -VR/32)Vu wh = Oh^R-^-W-V) +O l2 ( £ 500 ) 

By Lemma 3.3 and (6) we thus have 

(VR - VR/32)Vu = O Hr s(R^ 1+a - s ^)+O L ,(e 50a ). 

if T is chosen sufficiently large. By a Poincare inequality argument (adapted to 
the scale R), this implies that 

(Vr/2 Vr/w)(u -u) = 0^ +a -,(^- (1+a - 5)(1 - 5) ) + Hh (e 5m ), 

where 

_ _ Ib(0.R/2)\B(0,R/32) U 

U '- \B(0,R/2)\B(0,R/32)\ 
is the mean of u on the annulus B(0, R/2)\B(0, R/32). Since R = O(e" 20 °), we 
may replace the (e 500 ) error by O H i(e W0 ). Since u has bounded mass we 
see that u = 0(R~ 3 ^ 2 ), and it is then easy to check that {j] R / 2 — VR/ie)^ can be 
absorbed into the H i+«-s (R-( 1 + a - 5 )( 1 - s )) term. Thus we have 

(Vr/2 - VR/wh - HlR+a - S + 0h i (e 200 ). 

The contribution of the O h i (e 200 ) term to V((?7 fl /4 — 1!r/s)\ u \ 2u ) can easily be seen 
to be Ol 2 (^ 200 )) since u is bounded in this region thanks to (5). Thus by Cauchy- 
Schwarz the net contribution of this error term is 0(e 200 H/H2), which is acceptable, 
and so we may ignore this error and pretend in fact that 

On the other hand, from (5) we also have 

\\(Vr/2 ~ Vr/16)u\\l°° < -R -1 - 
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From an application of the fractional Leibnitz rule (sec e.g. [13], [25]) we thus have 

11(^/4 - Vr,b)\u\Mh?>-< Z R-tR-^-W-V, 
and the claim (66) follows by duality. ■ 

We can now verify (64) near the origin. This region has to be treated separately 
because the radial Sobolev bound (5) is too singular to be useful here, but the basic 
structure of the argument here will also be used later to prove (64) away from the 
origin. 

Proposition 8.3. If T is sufficiently large, then we have 

\(Vu wb (t)J) L ,\<\\f\\ H - a +e\\f\\ L 2 
for all test functions f supported on B(0, 1). 

Proof From the splitting (6) and Duhamel's formula (27) we have 

(Vu wb (i),/) L2 = (Vu wh (T>),e« T '-V A f) L2+l £ (VF(u(t')),e^-^ A f) L2 dt' . 

Letting T' — > +oo and using (8), it thus suffices to show that 

^ + °° \{VF(u(t')),e^'-^ A f) L ,\ dt' < \\f\\ H - a +e\\f\\ L i. (67) 

Let t := e~ 10 . We split this integral into the near future J t * +T and the distant 
future f£? T - To treat the distant future, we observe from (3), Leibnitz, Sobolev and 
Holder that 

\\VF(u(t'))\\ LHK ^<^ 

and hence by (19) we have 

{VF{u{t')),e^'-^f) L , < |t'-r 3/2 ||/|| Ll < |t'-t|- 3 / 2 ||/|| L2 

by the compact support of /. Thus this contribution is bounded by r -1 / 2 ||/||l 2 , 
which is acceptable by the choice of r. 

It then remains to control the near future J t * +T - We split F(u) into the portion 
rjiF(u) near the origin, and the portion (1 — rji)F{u) away from the origin. Let us 
first consider the portion away from the origin. By (65) with R = 1 we can bound 
this contribution to (67) by 

<| t+T ||^>-V(*'-*) A /|| jr -a +4 + e 100 ||e i (*'- t ) A /|U 2 dt'. 

The second term here is just T£ 100 ||/||i2 which is acceptable by choice of r. The 
first term can be controlled using (48) as 0(||/||^- Q ) (recall from the proof of 
Proposition 6.1 that when / is supported on B(0, 1) then we can take /+ = /) 
which is also acceptable. This concludes the treatment of the portion away from 
the origin. 
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It remains to show that 

f ' t+T \{V( m F(u(t'))),e^-^ A f) L2 \ dt' < \\f\\ H - a +e\\f\\ L i. 
t 

From the inductive hypothesis (64), with a replaced by a — 6, e replaced by e 100 , 
and R replaced by 1, we have by duality that 

Vu wh (t)=O Ha -s(l) + O L 2(e 10Q ) 

if T is large enough. From this and (7) we obtain in particular the local estimate 

m u wh (t) = O Hl+a -*(l) + O m (s W0 ) 

where i] 2 is a spherically symmetric cutoff to B(0, 1) which equals 1 on B(0, 1/2). 
From (6) and Lemma 3.3 we thus see that 

m u wh (t) = O m+a -*(l) + O H i(e W0 ) 

if T is chosen large enough. 

Let us use u er ror to denote the second term O^ife 100 ). This error term con- 
tributes a number of terms to V(r7iF(u)), which roughly look like either u 2 Vu error 
or uVuu error (there is also a term \u\ 2 u error arising from when the derivative hits 
the cutoff). In all cases we see from (3) and Sobolev that these expressions are all 
O L 6/5(e 10 °). On the other hand, by (20) e*(*'-*) A / = L * L% (||/|| 2 )- Thus by Holder 
the total contribution of the error is at most O(e 100 T 1 / 2 ), which is acceptable by 
choice of t. Thus we shall ignore the error term and pretend that 

\\ V2 u(t)\\ Hl+ <,-t < 1. (68) 

Let us now split up the integral J t * +T again, into the immediate future J t ' +1 and the 

medium-term future . For the medium-term future we use duality to rewrite 
this contribution as 

r +r |<e i(t -*' )A r?iF(w(i')),V/) i2 | dt'. 
Jt+i 

Because the kernel of e^ t_t ) A is smooth and decays like (t — i')~ 3 / 2 , it is easy to 
verify that 

for any K > 0. In particular we may control this expression in H a . Since 
|| |tt| 2 w(i')||i < 1 by (28), we thus see that this portion of the integral is bounded 
by WJWh-" as desired. 

Finally, we consider the immediate future. For this term we use the well-known 
Kato local smoothing estimate [38], [48] 

\\me^-^f\\ Ll H -. + i/. ([tim]xR .) < ll/ll ff - a( R3 } ; 

this allows us to estimate the contribution of the immediate future by 

< HV(^(«(0))ll^ BS -x/» ([ti4+1]xR . ) . (69) 

We may freely replace F(u) by F{rj2u). Heuristically, the fractional Leibnitz rule 
allows us to write 

V a - 1 /2 V (r ?1 |r ?2W | 2 r ?2U ) » r n (T l2 u) 2 X7 a+1 / 2 ( m u), 
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which heuristically allows us to estimate (69) by 

which is then acceptable by (31) and (68), if 5 is sufficiently small. This heuristic can 
then be justified by a rigorous application of the fractional Leibnitz rule, for instance 
by decomposing everything into Littlewood-Paley pieces, or by using the Coifman- 
Meyer theory of paraproducts; we omit the details. (Note that the presence of the 
L°f norm here is not dangerous, and in any event we have some surplus regularity 
to waste in the H a+1 / 2 component if one is concerned about cndpoints). ■ 

In light of the above proposition, it will suffice to verify (64) for R > 1. Accordingly, 
let us set R > 1, t > T, and let / be supported on R 3 \_B(0, R). Observe from 
Proposition 7.1 that 

\(Vu wh (t),f) L .\<R-i+ s \\f\\ L > 

and so we may assume 

1 < R < (70) 

since the claim (64) follows automatically otherwise from the above discussion. 
Indeed the above argument (combined with a smooth partition of unity) also shows 
that we may now assume that / is supported on the annular region {R < \x\ < 
£ - 1 /( 1 -< 5 )} j an d we will now do so. We also observe that it suffices to verify (64) 
when R is either equal to zero or an integer power of two, since the intermediate 
cases then follow automatically. In particular, by (70) there are only a finite number 
of values of R to consider. Thus in order to prove (64) it suffices to do so for a fixed 
value of R, since one can then let T be the supremum of all such times T obtained 
for individual values of R. 

We thus fix R, t, f and now apply Proposition 6.1 to decompose / = /_ + /+ + 
f smooth- The contribution of f amoo th can be estimated by (7), (44) as 

\{Vu wh {t),f smooth ) L ,\ < ||V/ smoot/l || L2(R3) < R-( 1+a ^-V\\f\\ k . a 

which is acceptable. 

We now consider the contribution of /_, arguing as in the proof of Proposition 7.1. 
From (45) we may assume that 

l(v e 4 * A u (o),/_)|< e ||./|| L2(R3) 

and 

l(V e 44A U+ ,/-)|< £ ||/|| i2(R 3 ) 

if T is chosen sufficiently large. Thus by (6) it suffices to show that 

l<V«(t),/_> ia - <V e ** A U (0),/_M < R-V+W-VMWb-. +e\\f\\ L2 . 

But by Duhamcl's formula (27) (as in the proof of Proposition 7.1) we may estimate 
the left-hand side here by 

< /Vw')),e <(t '- t)A /-M dt'. 
Jo 
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Let t := e 10 . We shall split this time integral J Q into the distant past J Q T 
and the recent past J* . Note that i — r is positive if T is chosen sufficiently 
large depending on e. The contribution of the distant past is treated similarly 
to the corresponding contribution of the distant future in (67); we pick up an 
additional factor of 0(e _3 / 2 ( 1- ' 5 )) because / is now supported on the ball of radius 
Oie^ 1 ^ 1 ^^) instead of O(l), but this still does not affect the argument significantly 
because r is so large. 

It remains to consider the recent past. We have to show that 



/ |(V^( W (t')), e i(f '-* )A /-) i2 1 dt' < .R-d+^C 1 -^ ll/H + e ||/|| i2 
J t-r (71) 



Let r?_R/8 be a smooth cutoff to B(0,R/8) which equals 1 on B(Q, R/16). We split 
F(u(t')) into the part rj R /g,F(u(t')) near the origin and the part (1 — rj R / 8 )F{u(t')) 
away from the origin. Let us first deal with the part near the origin. For this part 
we use Cauchy-Schwarz, observing from (3) and Sobolev that F{u(t')) is in L 2 (R 3 ), 
to bound this portion by 



|Ve i (*'- t ) A /_|| L2(s(0)il/8) dt'. 



But this component is acceptable by the exponential decay bound (47) near the 
origin (choosing a, (i in (47) sufficiently large). 

It remains to deal with contribution of (1 — r] R / s )F(u(t r j) to (71). By Lemma 8.2 
we see that this contribution 

< R - { l +a -W- S) ^ \\ {{x) -2 e i(t'-t)A f _ llH _ a+s +e 100|| e (t'-t)A / _|| ia) dt , 

J t — T 

The contribution of the e 100 term is at most 



< £ 10( V||/- 



is- 1 ' 



which is acceptable by (43) and the definition of r. Now consider the main term. 
Applying (49) we may bound this by 

< J j-(l + a-5)(l-«) J j-l+«|| / _|| ff _ a 

which is acceptable. This proves (71) and concludes the treatment of /_. The 
treatment of /+ is very similar, and involves evolving forward in time instead of 
backward in time (cf. the proof of Proposition 8.3), and so uses u + instead of u(0) 
but is otherwise almost identical. (Note also that the integration parameter t' is 
now greater than t and hence also greater than T, so there is no difficulty ensuring 
that t' is large). This proves (64), which closes the induction and completes the 
proof of the Theorem. ■ 



We now wrap up by using (64) to prove Theorem 1.3. We first undo the duality 
in (64). Let a > 0, S > and e > be chosen later, and suppose that t is 
sufficiently large depending on a, 5, e; our implicit constants can depend on a, 5. 
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For the portion of u w ^ near the origin, we observe from (64) (replacing e by e 100 ) 
that 

\( V2 Vu wh (t)J) L i\<\\f\\ H - a +e W0 \\f\\Li 

for all test functions / and for any a > 0, e > 0. By duality, this means that we 
may decompose 

V2Vu wh (t) = O ffQ(R3) (l) + L2{K3) (e im ). 

On the other hand, from (7) we have that r]2U w 'fo{t) — O h1 ^3^(1). A Poincare 
inequality argument (writing 'i]i/2 u w \ > (t) as a localized fractional integration oper- 
ator of r^Vu^W of order -1, plus an infinitely smoothing operator of T]2U w \)(t)) 
then easily yields that 

»h/2ttwb(*) = °h^(R 3 )( 1 ) + °m (R 3 ) ( £l °°)- 

In particular, we can find a function u^ (t) for all sufficiently large times t such 
that 

||^ w bW-%,ollHi(R 3 )^ £l0 ° ( 72 ) 

and 

H u b,oll ff c +1( R^ < 1- 

By applying r\i if necessary we may ensure that Q is supported on 5(0, 2). From 
Sobolev embedding we see that 

|V J 'u bi0 (x)| < 1 (73) 

for all < j < a + 1/2. 

This is enough to construct near the origin. We now start working in dyadic 
shells outside of the origin. We first observe from Proposition 7. 1 that it will suffice 
to construct on, say, the ball B(0, £ -10 ), since the energy of u w ^ outside this 
ball is asymptotically zero. Accordingly, we let 1/4 < R < e~ 10 be a power of 2. 
From (64) we have that 

|<(»?3 2 * - »7 2 «)Vu wb (t),/> ia | < iz-^x 1 -^!!/!!^-- +e 100 ||/IU=. 

for all test functions /, if t is large enough (note that there are only a finite number 
of R involved, so we can use a single time threshold t >T = T(e, u, a, S) for all of 
them). By duality again, this implies 

(V32R ~ V2R)Vu wh = 0^(i?- (1+ " )(1 - 5) ) +O L2 (£ 10 °). (74) 

We now claim that a Poincare inequality argument gives 

(Vsr - V4 R )u wh = O^iR-^+^-^ + O^ is 100 ). (75) 

This time we shall perform this argument in detail, taking note of the scaling factor 
of R. Wc begin by observing from integration by parts that for any function / on 
R 3 and any unit vector u, we have 

/>OC 

f(x) =- T] R (r)bj ■ Vf{x + ruj) + T]' R (r)f(x + ru>) dr. 
Jo 
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Averaging this over all u> G S 2 and then removing the polar co-ordinates, we obtain 
m = C | r3 V(y)j^ ■ V/(x + y) + ^f(x + V) dy. 

Applying this to / = u w ^ when x lies in the support of f] S R — r]4R, we see in 
particular that 

(?78K - ViR)u wh (x) =C(r] 8R (x) - r} 4R (x)) ^ 3 rj R (y) ■ V((r) 32 R - V2R)u wb ){x + v) d V 

+ C(t) S r(x) - t]4r{x)) r}R(y)\y\~ 2 {{V32R - r) 2R )u wh )(x + y) dy. 
jR* 

We now substitute (74) into the first term, and (7) into the second term. Observe 
that convolution with r}n{y)-A^ is a smoothing operator of order -1 and maps H R 

to H R +1 with constant O(l) (the factor of R can be easily scaled out) and L 2 to 
H R , again with a bound of O(l). Meanwhile, convolution with T] R {y)\y\~ 2 is a 
smoothing operator which maps L 2 to H R +1 with constant 0(1). The claim (75) 
follows. 

From (75) we may find a function R (t) for each 1/4 < R < e~ 10 and all suffi- 
ciently large times t such that 

\\(V8R - V4 R )u wh - u^rWh^h^ < £ 10 ° (76) 

and 

By applying t)iq R — r\2R to R if necessary we may assume that R is supported 
in the annulus B(0, 16R)\B(0, R). From Sobolev embedding (rescaled by R) we 
thus see that 

\V j u hiR (t,x)\<R-l-jR^ 5 (77) 
for all integers < j < a + 1/2. If we now define 22 

u h ~ u b,o + 2 u b,R 

l/4<fl<e- 10 

where R ranges over powers of 2, we thus see by adding up (73), (77) that 
obeys the symbol bounds (14), if a is chosen sufficiently large depending on J, and 
S chosen sufficiently small depending on a. Furthermore, by telescoping (72) and 
(76) we see that 

ll%flo«wb(*) _ u b(*)Hiji(R 3 ) ~ £ 10 °logl/e, 

where i?o is the largest power of two less than or equal to e~ 10 . But from Proposition 
7.1 and Hardy's inequality we see that 

||(l-»»flo)«wb(*)lljji(R 8 ) ~ £ 



22 Strictly speaking, this only defines for sufficiently large times. But for bounded times we 
can set arbitrarily to (say) 0. In practice we will actually need to take a sequence of values of 
e tending to zero, with corresponding time thresholds T e > tending to infinity, and redefine 
according to a new value of e every time we cross one of the time thresholds, in order to make the 
error in (13) o J j 1 (l) rather than merely Of I1 {e) for any given s. We omit the details. 
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if t is sufficiently large. Thus we obtain 

u wb(*) = %(*) +°Hi( £ ) 
which will imply (13) since £ can be arbitrarily small (though of course this forces 
t to be increasingly large; see previous footnote). Finally, from (76) we see that 

\\(V8R ~ V4ii)u wh - % !i? ll H1( R3) < #£ 10 ° 

which telescopes (together with (72)) to obtain that 

ll'J8flo«wb - %Hhi(R 3 ) ~ ^oe 100 - 

Since Rq < e -10 , the claim (43) then follows from (7). Finally, we can easily ensure 
that itu is spherically symmetric (either by making sure each step of the argument 
preserves spherical symmetry, or else averaging over rotations at the end of the 
argument. This concludes the proof of Theorem 1.3. ■ 



9. A DIGRESSION ON VIRIAL-TYPE IDENTITIES 



Before we begin the proof of Theorem 1.4, we pause to review the derivation of 
virial identities such as (16). First observe that if u is any solution to (1), then we 
have the momentum flux identity 

d t Im( Uj u) = -2d k Rc(uJu k ) + idjA(M 2 ) + ^>| 4 (78) 

where j, k range over 1,2,3 with Uj := djU = d Xj u and the usual summation 
conventions; this can be verified by direct calculation (or see [45]). Integrating this 
against a suitable gradient ctj = dja for some explicit weight a(x) to be chosen 
later, we (formally) obtain after some integration by parts 

d t J ajlm(ujH) dx = J —2ajdkRc(uJuk) + ^ajdjA(\u\ 2 ) + ^ajdj\u\ 4 dx 

2a jk Re(uJu k ) - -(AAa)|w| 2 - -(Aa)|u| 4 dx. 

Thus for instance if a(x) := \x\ 2 , then we obtain Glassey's virial identity (16). 
Unfortunately we do not have quite enough decay on to make the left-hand 
side of (16) finite, and so we will instead use a slightly different weight a, namely 
a(x) := \x\ 2 (x)~ s rm, where S > is a small exponent and R ^ 1 is a large radius, 
and rj R is a cutoff to B(0, R) which equals 1 on B(0, R/2). Observe that 

AAa = 0{S)(x)- 2 - s + 0(R- 2 - s ) 

and 

Aa = (6 + 0(S)){x)- s + 0(R- S ) 
on the ball B(0,R), while 

a jk Re(u]u k ) = ((1 + 0(5))(x)- 5 + 0(R- s ))\Vu\ 2 . 

We thus have 

d t I 2ajlm{ Uj u) dx = [ 4:(x)- s \Vu\ 2 -3(x)- 5 \u\ 4 dx+0((S+R- s )\\u\\ 2 H1 + \\u\\ A Li ). 
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We of course can estimate ||w||l 4 by ||u||#i by Sobolev embedding. Integrating this 
in time, and using the fact that |Va| = 0((x) 1_<5 , we obtain 

|- [ T+T [ A(x)- s \Vu\ 2 - 3<x)- 5 M 4 dxdt\ < 
T Jt Jb(q,r) 

- / (xf-'MQWVuWl dx 

T t=T,T+r J B(0,R) (,' y j 

+ (S + R- S )(\\u\\ LrHimT+T]xR s ) )\ 
This will be our main tool in proving Theorem 1.4, which we do next. 



10. Proof of Theorem 1.4 



We now prove Theorem 1.4. Let u be as in the theorem, and let e be arbitrary. We 
need to show that for t > sufficiently large depending on e, and T > sufficiently 
large depending on e, r, that 

l r It +T Ik 3 4|V%|2 ~ 3|%|4 dxdtl ~ £ ' (80) 



Fix e. We now choose a number of parameters; it will be important to pay attention 
to the order in which these parameters are selected. For reasons which will be 
apparent later we will need a large radius R > depending on e > to be 
chosen later. To apply (79) we shall need a small exponent 5 > depending on 
e, i?o to be chosen later, as well as an even larger radius R » Ro depending on 
S, e, Rq to be chosen later, and then we will choose r sufficiently large depending 
on e, <5, i?o, R- Finally we let T be a sufficiently large time (depending on u and all 
previous parameters) again to be chosen later. Let m w k t be the solution to (1) 
with initial data u ^ T (T) = u w ^(T). Similarly define u^ T . We now show that 
these solutions will approximate the actual functions on the time interval 

[T, T + t] if T is large enough. 

Lemma 10.1. IfT is sufficiently large (depending on t), then we have 
\\u wb (t)-u wbT (t)\\ Hl <s w R- w 

and 

h b (t)-u bT (t)\\ H1 <s w R- w (81) 

for all t € [T,T + t]. In particular, from the H 1 local well-posedness theory we see 
that u bT and u wbT are well-defined for all t G [T,T + r] (note that the L 2 norm 
ofu bT stays bounded by mass conservation). 



Proof Fix T, r; all spacetime norms will be on the slab [T, T + t] x R 3 . We will 
assume a priori that ml and m w k exist and are in H 1 on the entire time interval 
[T, T + t], this a priori assumption can then be removed by the usual continuity 
argument, e.g. letting r increase continuously from zero. 
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We begin with u w ^. Write v = u w ^ — u w ^ T1 then v(T) — and we have the 
equation 

iv t + Av = F(u wh ) - F(u wh -v) + {id t + A)u wb - F(u wh ). 

However, from (10) we have that 

\\(id t + A)u wh -F(u wh )\\ L i H i <e' 

for an e' > (much smaller than e) which we will choose later, if T is sufficiently 
large depending on e', r, u. We can thus write the above equation as 

iv t + Av = F(u wh ) - F(u wh -v) + L i H i (e'). 

Let X denote the norm 

\\v\\x ■= \\v\\l°°hi + IMIl^l-- 

From the Strichartz inequality (24) (and energy estimates and Duhamel's formula) 
we have 

Hx < \\v(T)\\hi + II^Kb) - ^Kb - v )hlHi + e'; 
we will refrain for the moment from using the fact that v(T) = 0. However, from 
Leibnitz and Holder one easily checks that 

\\F(u wh ) - F(u wh - v)\\ L , Hi < \\v\\ x (\\v\\x + IKblU) 2 

(basically, one applies the derivative in the H 1 norm to the cubic expression F(u w y ) ) — 
F(u w ^ — v), and whatever the derivative hits is placed in Lfl? xl and the other two 
factors in LfL"^). Thus by Holder we have 

HI* < HT)\\ Hk + r^\\v\\ x (\\v\\ x + \\u wh \\ x f + e>. 

Now observe that Lemma 3.2 gives bounds on the X norm of u, which are also 
obeyed by the free solution e ltA u + . Hence by (6) we have Hu^Hx < (t) 1 / 4 , thus 

\\v\\x < \\v(T)\\ hi +t^ 2 \\v\\ x (\\v\\ x + (r) 1 / 4 ) 2 +s'. 

This estimate does not help us directly when r is large. However, if we decompose 
the interval [T, T + r] into 0(t) intervals [ti, t\ + r'] of small length t' « 1, and 
apply this estimate on each interval, we may show (using continuity arguments in 
the usual manner) a bound of the form 

ll W llx([t lltl+ r']xR 8 )^IK*l)ll^+ e ' 

assuming that was sufficiently small. In particular we can control the 

H 1 norm of v(ti + t') in terms of that of v(t\), plus an error of e' . Iterating this 
we see that (if e' is sufficiently small depending on r', e, R) we have 

sup \\v(t)\\ m <e w R- w 
te[T.T+r] 

as desired. 

Now we handle u^. From (13) and the bounds just proven, it will suffice to show 
that 

sup \\u hT (t)-u whT (t)\\ m <e w R- w 
te[T,T+ T ] 
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for T sufficiently large. Write w := u bT (i) — u wbT (f). By (13) we may ensure 
that 

\\<T)\\fr <e' 

for any e' > to be chosen later, if T is sufficiently large depending on e' , u. Since 
i* b j, and w wb T both obey (1), we see that w obeys the equation 

iw t + Aw = F(u wh T + w)- F(u wh T ). 

We now perform a homogeneous version of the previous analysis. Let X denote the 
norm 

Mllx := MILfffi + MIIl|£~j 

then by Strichartz again (in particular (24)) we have 

Mix £ IH T )llffi + ll F Owb,T + w ) - ^("wb^lliJiTi- 
By Leibnitz and Holder we have 

\\F(u wh T +w)- F(u wh T )\\ Ll ^ < T 1/2 |Mlx(|| UwbiT ||;f + \Mx) 2 - 

But from the previous discussion we already have bounds of the form ||u wb || x ~ 1 
when r is small enough. Thus if we subdivide [T, T + r] into small intervals and 
choose e' sufficiently small as before, we may obtain 



sup Hljj, <e 10 i?- 10 

*G[T,T+r] 



as desired. 



We now apply (79) to u b T to obtain 
1 f T+T f 

|- / / 4(x)- ,5 |Vw bT | 2 -3(a;)- d |u bT | 4 dxdt| 
r Jt Jb(o,r) 

<- sup / (a;) 1 - 5 | %T (i)||V UbT (t)| 

=T,T+t Jb(0,R) ' ' (o^J 



T t= 

-<5\ 



+ ( ( 5 + i?- d )(|| >T || Lr ^ ([TT+r]xR3) ) 4 . 
We now estimate the various terms in this expression. 

First of all, from (15) we have 

\\u hT (T)\\ m = \\u h (T)\\ m <l 
and in particular by mass conservation 

ll M b,TlL~L2([T,T+r]xR 3 ) ~ L 

Also from (81) we have 

H M b,TllL~Hi([T,T+r]xR 3 ) ~ 1 - 

Thus we have 

(H M kTllL~iJi([T : T+T]xR 3 )) ~ 1 

and thus (if we choose 5 small enough depending on e, and i? large enough depend- 
ing on 5) 

(5 + ir 5 )(|| >T || Lr ^ ([TT+r]xR3) ) 4 < e. 
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We now consider the first term in the right-hand side of (82). Crudely bounding 
(x) 1 - 5 by R and using the above H 1 bounds on u b T we see that this term is 
bounded by 0(R/t). Thus if we choose r large enough depending on R then this 
term is also 0(e). We have thus obtained 

i r T+T r 

\- / 4(x)- <5 |Vu bT | 2 -3(x)- s \u hT \ 4 dxdt\ <e. (83) 

T JT JB(0,R) 

We now divide the spatial region of integration B(0,R) on the left into the part 
B(0,Rq) near the origin, and the part B(Q, R)\B(0, Rq) away from the origin. 
First consider the part away from the origin. From the bounds in (14) (or from 
Proposition 7.1), applied for instance with 5 = 0.1, we have that 

\[ |V % | 2 + | % | 4 dx|<i? 2+0 - 2 (84) 

JR \B(0,Ro) 

if T is sufficiently large. Using (81) and Sobolcv and Holder, we thus have 

|v %T | 2 + K T | 4 d*|<i? - 2+0 - 2 

iB(o,n;)\B(o,Ro) 

if T is sufficiently large. In particular, if Ro is chosen sufficiently large depending 
on e we have 



Jbi 



Ib 



|Vw bT | 2 + |u bT | 4 dx < e. (85) 

I B(0,R)\B(Q,R a ) 

On the other hand, inside the ball B(0,Rq) we have (x)~ s = 1 + O(SR ) (for 
instance), and so again by the H 1 bounds on « b T we have 

f 4(x)^|Vu bT | 2 -3(x)-' 5 |u bT | 4 dx= f 4|Vii bT | 2 -3|u bT | 4 dx+(9((5i?o). 

JB(0,Ro) ' JB(0,Ro) 

Again, if we choose 5 sufficiently small depending on R then 0(5 Ro) can be made 
less than e. Combining this with (84), (85) we thus have 

Mx)- 5 \Vu h T \ 2 - 3(x)- 5 \u h T | 4 dx = f 4|V % | 2 - 3| % | 4 dx + 0(e). 
Combining this with (83) we obtain (80) as desired. ■ 



Remark. The above argument can be refined slightly to show that 

lim sup | / / 4|Vu b | 2 -3|u b | 4 dxdt\ < (t) ct 

for any a > 0, with the implicit constant depending on tr. This can be done by 
setting Rq, R to be large powers of t, and S to be a small power of r, and using 
(14) to improve the crude bound of 0(R/t) obtained above for the first term of 
(82) to O{R 0+ /t). Wc omit the details. 

Remark. By choosing different multipliers a (not just perturbations of \x\ 2 ), it is 
likely one can prove other asymptotic identities of this type. For instance, it seems 
likely that one has the asymptotic Morawetz identity 

i r T+T r \u u \ 4 

lim lim sup |- / (47i> b (i,0)| 2 - / dx)\ = 0, 

r^+oo T^+oc T J T u JR 3 \X\ 
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which (formally at least) is associated to the multiplier a = \x\, although the singu- 
larity at has to be treated here with some care. More generally, it seems that the 
momentum flux (78), when applied to uu, should converge asymptotically to zero 
in some weak sense when averaged over increasingly large intervals of time. If one 
could get stronger control on the decay of the momentum flux, e.g. if it converged 
in a suitably strong sense to zero without any averaging in time, this would be 
substantial progress towards the soliton resolution conjecture (solitions have zero 
momentum current, and are probably the only spherically symmetric solutions to 
(1) with this property). However, to do this it seems one would first have to dis- 
cover a technique that would eliminate the possibility of breather solutions (smooth 
localized periodic solutions to (1) which are not of the form u{t,x) = Q(x)e lwt ). 
We do not know of a way to rule out such solutions; the monotonicity formulae aris- 
ing from virial-type identities do not appear appropriate for this task. Note that 
the completely integrable ID equation supports some breather solutions (see [52]), 
but this may be a phenomenon caused by complete integrability and we do not 
know if there is a similar phenomenon in higher dimensions, say in the spherically 
symmetric case. 
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